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1.  INTRODUCTION 

1.1  General  Objectives  and  Scope  of  the  Program 

The  general  objectives  of  the  program  reported  in  this 
document  were  the  implementation  and  support  of  research 
directed  toward  measurements  of  the  energetic  charged 
particle  environment  in  near-earth  space. 

Specific  objectives  were:  the  design,  construction  -.and 

calibration  of  instrument  packages  for  the  U.S,  Air  Force 
Satellites  DMSP  (Flights  #2,  3,  4  and  5)  and  P78-1,  and  field 
support  for  the  same  missions;  the  study  of  the  response, 
through  computer  and  laboratory  simulations,  of  existing 
proton-telescopes  to  the  natural  radiation  environment;  the 
data  analysis  of  proton  data  obtained  from  previous  satellite 
missions;  finally^he  preliminary  study  and  design  of  a 
charged  particle  telescope. 

1.2  Organization  of  the  Report 

Section  2  describes  the  instruments  which  were  flown 
on  the  DMSP  and  P78-1  satellites.  This  section  also  con¬ 
tains  a  brief  description  of  the  objectives  of  each  mission, 
and  operation  summary.  Some  sample  data  are  also  presented. 
Photographs;  block  diagrams  and  first  launch  engineering 
evaluations  are  included  where  applicable. 

Section  3  describes  the  simulations  made  to  model  the 
reponse  of  the  proton  telescopes,  LEPS,  to  the  proton  en¬ 
vironment.  These  proton  telescopes  were  flown  on  U.S.  Air 
Force  satellites  S3-2  and  S3 -  3. 
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In  Section  4,  the  data  analysis  activities  are  presented. 
These  include  the  data  handling  software  development,  some 
results  of  the  analysis  of  data  from  the  S3-2  satellite,  and 
the  evaluation  of  the  generalized  geometric  factor  of  a 
particle  telescope. 

Section  5  contains  a  summary  of  the  efforts  toward  the 
design  of  a  proton  spectrometer  with  a  high  sensitivity. 

Section  6  contains  a  report  (specifically)  of  research 
in  the  detection  of  protons  in  CR-39  plastic,  aimed  at 
interest  in  the  Long  Term  Duration  Mission. 

In  Section  7  various  facets  of  software  development 
for  the  MAGAF  system  are  discussed  and  some  programs 
are  listed. 

The  use  of  data  from  the  MAGAF  network  in  the  light 
of  micropulsation  studies  is  discussed  in  Section  8. 
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2. 


INSTRUMENTATION 


2.1  The  Defense  Meteorological  Satellite  Program  (DMSP) 

The  Defense  Meteorological  Satellite  Program  (DMSP)  con¬ 
sists  of  a  series  of  satellites  in  450  nautical  mile,  polar 
orbits,  each  carrying  a  payload  of  meteorological  sensors. 

In  addition  to  the  primary  payloads,  there  are  on  each  of 
the  satellites,  several  instruments  one  of  which  is  a  small, 
light  weight  (3  lb.)  electrostatic  analyzer,  the  SSJ/3. 

2.1.1  The  SSJ/3  Payloads 

The  purpose  of  these  payloads  is  to  measure  precipitating 
electron  fluxes  over  the  polar  caps,  in  the  auroral  zones  and 
at  the  plasma  pause  regions.  These  data,  in  particular  the 
data  acquired  in  the  auroral  zones  during  magnetic  storm  and 
substorm  periods,  will  be  analyzed  together  with  data  obtained 
by  other  sensors  such  as  optical  sensors  on  board  the  DMSP 
satellites,  and  the  AFGL  magnetometer  networks.  The  outcome 
of  such  an  analysis  will  be  an  improvement  upon  existing 
substorm  models,  the  understanding  of  the  mechanisms  and 
sources  of  particle  acceleration,  and  the  possibility  of  pre¬ 
dicting  magnetic  activity. 

Each  SSJ/3  payload  consists  of  a  set  of  two  electrostatic 
analyzers  which  perform  a  differential  energy  analysis  of  low 
energy  electrons  in  the  energy  range  50  eV  to  20  keV.  The 
analyzers  are  of  the  curved  plate  type  and  utilize  channel 
electron  multipliers  to  detect  electrons  of  selected  energy. 

A  16-point  spectrum  is  produced  each  second.  A  sun  sensor 
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turns  the  instrument  off  when  the  sun  is  in  the  field  of  view 


prolonging  the  life  of  the  channeltron  electron  multipliers. 

A  detailed  description  of  the  SSJ/3  has  been  published 
elsewhere  (AFG L- TR- 77- 0  1  20)  .  Table  2.1  lists  a  summary 
of  the  instrument  characteristics. 

Two  independent  calibrations  have  been  performed  on  the 
SSJ/3  packages.  First,  before  the  first  unit  had  flown,  a 
calibration  was  done  at  the  Aerospace  Corporation.  A  com¬ 
parison  was  made  between  this  calibration  and  a  Monte  Carlo 
computer  simulation  of  the  instrument.  The  two  were  in 
excellent  agreement.  The  results  of  this  calibration 
have  also  been  published  ( AFGL- TR- 7 7- 0 1 0 2 )  . 

Because  of  the  wide  interest  in  the  data  from  these 
instruments,  a  second  calibration  was  performed  recently  on 
units  4  and  5  using  a  more  sophisticated  beam  facility  at 
Rice  University,  The  results,  to  be  published,  confirmed 
the  previous  results  and* also  provided  data  on  UV  sensitivity 

2.1.2  Summary  of  Operations 

Four  SSJ/3  packages  were  provided  for  DMSP  block  5D 
flights  #2,  3,  4  and  5. 
a)  Flight  #2 

This  satellite  was  launched  June  1977,  and  our  instru¬ 
ment  has  been  operating  properly  ever  since.  Data  from  this 
package  has  been  reliable  and  useful  to  the  DMSP  effort. 
Typical  data  is  shown  in  Figure  2.1  and  2.2. 
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Table  2  .  1 


Particles  Detected 
Detectors  Used 

Method  of  Energy  Analysis 
Number  of  Energy  Bins 

Energy  Bins 


Electrostatic  Analyzer 
SSJ/3 

Summa  ry  of  Characteristics 
E 1 ectrons 


2  Channeltron  Electron  Multipliers  for 
each  of  the  two  sets  of  ESA  plates 

Voltage  Stepping  on  ESA  plates 

8  for  each  set  of  ESA  plates; 

Total  of  16 


Large  1 

Plates 

Sma  1 1 

Plates 

20,000 

eV  +  3% 

1  ,045 

eV  +  32 

13,700 

+  3% 

661 

+  3% 

8,990 

+  3% 

434 

+  3% 

5,500 

+  3% 

264 

+  3% 

3,790 

+  3% 

183 

+  3% 

2,290 

+  3% 

no 

+  3% 

1  ,590 

+  5% 

77 

+  5% 

1  ,060 

+  6% 

51 

+  62 

Acceptance  Angles 

Large 

ESA: 

1.6°  FWHM  across  the  apertures 
8,0°  FWHM  along  the  apertures 

Sma  1  1 

ESA: 

3.7°  FWHM  across  the  apertures 
4.8°  FWHM  along  the  apertures 

Normalization  Constants 

Large 

ESA  : 

-  4  2 

1.3x10  cm  -ster. 

AE 

Sma  1  1 

ESA: 

4. 3xl0_^cm2-ster. 

E 

Large 

ESA  : 

4.0? 

Small 

ESA: 

7.22 

Data  Rate 

One 

omplete  Spectrum  per  second 

Table  2,1  (cont  . ) 


Dwell  Time  at  Each  98  ms. 

Energy  Level 

Digital  Data  Format  (9  bits  per  channel)  x  (1G  channels) 

=  144  Bits 

Note :  Each  9  bits  in  every  channel  are  in  logarithmic  form,  the 
five  least  significant  bits  being  the  mantissa  and  the 
remaining  four  the  exponent.  This  number  is  converted  to 
decimal  form  according  to  the  following  relationship: 


N  =  2y(x+32)  -  33 


Where;  ,0  0  0  0,,0  0  0  0  0 

V 

y  x 


^LSB 


Data  Readout 


y  ^  ^ 

^ooo^boooq) 

51  e  V 


The  first  bit  to  be  read  out  is  the  least 
significant  bit  ( L  S  B )  of  the  highest 
channel  followed  by  the  next  to  the  highest 
one  etc.  ,  i . e. : 


.  'O  0  0  00  0000,0  00000  oocr, 

13,700  eV  20,000  eV 


1st  bit  out 


[ . 144  bits . ] 

read  out  in  one  group 

at  the  end  of  each  second 


Analog  Monitors 


Size 


Plate  voltages:  5.0  volts  to  .  25v 
Power  supply  :  2.5v 
Temperature  :  2.5v  at  room  temp. 

5 . 50 i n  X  3 . 29 i n  X  5.  lOin 


Weight  3. 04 Gibs 

Power  Dissipation  .125  Watts 


o 


POLEWARD  POLAR  CAP 

BOUNDARY  ARC 


JI  1  /CKSIC  i 


SSJ/3 


I 


16-POINT  ENERGY  SPECTRUM  OVER  THE  POLAR  REGION  FOR  10  CONSECUTIVE  SECONDS 
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b)  Flight  #3 

Launch  date  was  June  1978.  Our  instrument  on  this 
satellite  functioned  for  a  few  months,  but  an  I.C.  failure  in 
the  regulator  circuit  forced  us  to  shut  down  the  instrument 
when  the  H,V.  stepping  circuit  became  inoperative. 

c)  Flight  #4 

Launch  date  was  June  1979.  We  just  received  the  first 
"quick  look"  data  from  this  spacecraft  indicating  that  our 
instrument  is  working  properly. 

d)  Flight  #5 

Classified  launch  date. 

The  data  recorded  by  these  instruments  has  been  recog¬ 
nized  as  being  of  extremely  high  quality,  free  from  noise  and 
other  spurious  effects  which  have  plagued  previous  instru¬ 
ments  of  the  same  general  type  (Fig,  2,1,  2,2).  Because  of 
its  exceptional  quality  as  well  as  the  great  current  interest 
in  precipitating  electrons  on  the  part  of  the  geophysics 
community,  the  data  are  being  used  by  a  number  of  agencies 
and  research  groups,  including  AFGL,  Applied  Physics  Labora¬ 
tory  of  John  Hopkins,  NOAA,  and  the  Aerospace  Corp,  which  is 
correlating  our  data  with  their  X-ray  data,  in  addition  to 
the  Air  Force  Air  Weather  Service  which  is  the  primary  user, 

2.2  The  P78-1  Satellite 

The  main  payload  on  the  P78-1  satellite  consists  of  the 
Advanced  Research  Projects  Agency  (ARPA) ,  301  Gamma  Ray 
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Spectrometer.  There  are  six  secondary  payloads,  one  of  which 
is  CRL-251,  High  Altitude  Particle  Spectrometer. 

2.2.1  The  CRL-251  Payload 

The  mission  of  this  payload  is  to  measure  the  flux  of 
particles  precipitating  over  the  polar  cap  regions,  the 
auroral  zones,  and  in  the  plasma  pause,  and  to  provide  boun¬ 
dary  conditions  for  the  particle  population  of  the  magneto¬ 
sphere.  These  measurements  provide  data  on  the  link  between 
observed  trapped  and  interplanetary  populations,  and  the 
particles  entering  the  atmosphere. 

In  addition,  the  payload  will  acquire  electron  data  in 
the  auroral  zones,  particularly  during  magnetic  storm  and 
substorm  periods.  The  results  of  these  studies,  together 
with  existing  substorm  models  and  with  data  from  other 
sources,  will  be  used  to  improve  upon  existing  substorm 
models.  This  will  include  f i e  1  d -a  1  i gn ed  currents  and  the 
mechanisms  and  sources  of  particle  acceleration.  Further¬ 
more,  these  data  will  be  correlated  with  the  AFGL  magnetic 
network  data  to  help  AFGL  in  their  efforts  toward  predicting 
magnetic  activity.  The  results  should  be  useful  to  the  ex¬ 
periments  investigating  high  energy  particles  and  the  auroral 
regions  including  the  primary  mission  objectives. 

The  CRL-251  payload  consists  of  a  set  of  four  electro¬ 
static  analyzers  wnich  perform  a  differential  energy  analysis 
of  low  energy  electrons,  simultaneously  in  two  look  direc¬ 
tions.  A  16-point  spectrum  is  produced  for  each  look 
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direction  every  256  msec  by  the  low  and  high  energy  analyzers. 
The  analyzers  are  a  conventional  curved  (cylindrical)  plate 
design  with  balanced  positive  and  negative  deflection  volt¬ 
ages  . 

The  CRL-251  sensor  box  is  located  in  Compartment  9  of 
the  wheel  assembly  (Figure  2.3).  The  box  is  oriented  so  that 
the  look  axis  for  one  set  of  analyzers  is  2.5  deg  from  the  +X 
axis  of  the  wheel  and  the  look  axis  for  the  other  set  of 
analyzers  is  2.5  deg  from  the  -Y  axis  of  the  wheel.  The  spin 
axis  is  perpendicular  to  the  orbital  plane  so  that  the  wheel 
rotates  in  the  plane.  Thus,  the  sensors  look  toward  and 
away  from  the  earth  during  each  spin.  Since  the  particle 
pitch  angle  is  the  angle  between  the  particle  direction  and 
the  magnetic  field  vector,  the  sensors  go  through  360  deg. 
of  pitch  angles  for  the  particles  collected,  with  a  good 
angular  resolution. 

This  instrument  was  designed  to  collect  electrons  in  the 
energy  range  of  50  eV  to  20  keV.  Two  sets  of  curved  parallel 
plates  have  a  time-sequenced  variable  electrostatic  field  to 
deflect  the  electrons  toward  channeltron  multipliers  which 
are  used  as  detectors.  One  set  of  plates  collects  electrons 
in  the  energy  range  of  50  eV  to  1  keV,  while  the  second  one 
collects  electrons  from  1  to  20  keV.  These  are  swept  simul¬ 
taneously  to  provide  four  complete  spectra  every  second, 

F.ach  set  of  plates  is  divided  into  eight  energy  bins,  so  that 
each  set  of  spectra  consists  of  16  energy  bins. 


THIS  PAGl?  IS  BPS  I  QUALITY  TRACTI  CARLS 
FKIUM  (A>i  i  x  £0  IXL#C 
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digital  solar  aspect  sensor  elec 

DIGITAL  SOLAR  ASPECT  SENSOR 


2.2.2  Operational  Description 

Figure  2.4  shows  the  payload  sensors  with  the  cover  re¬ 
moved.  The  large  slits  are  the  high  energy  analyzers  (1  to 
20  keV)  for  each  set.  The  middle  slits  are  for  the  sun  sensors 
and  the  other  slits  are  for  the  low  energy  analyzers  (50  ev  to 
1  keV) .  The  apertures  on  the  case  are  directly  in  front  of 
these  slits  or  particle  collimators  (Figure  2.5)  and  behind 
the  collimators  are  the  curved  cylindrical  deflection  plates. 
Different  voltages  are  applied  across  these  plates  in  se¬ 
quence  to  allow  electrons  of  different  energies  to  pass 
through.  The  plate  separation  and  the  applied  voltage  deter¬ 
mine  the  energy  band  of  the  collected  electrons  such  that 
the  mean  energy  is  directly  proportional  to  the  applied  vol¬ 
tage  difference  and  inversely  proportional  to  the  plate  separ¬ 
ation.  The  sector  angle  of  the  plates  and  their  separations, 
together  with  aperture  and  exit  slits,  determines  the  counting 
efficiency  of  the  sensors. 

Additional  baffles  in  the  aperture  assembly  eliminate 
sunlight  (UV)  scattered  in  the  detectors. 

Two  channeltrons  are  mounted  behind  the  exit  aperature 
in  each  sensor.  The  electrons  impinge  on  the  cone-shaped  area 
and  produce  secondaries  which  cascade  through  the  channeltron 
producing  a  detectable  pulse  at  the  end  which  is  then  processed 
by  the  internal  electronics.  The  operating  lifetime  of  the 
channeltrons  is  extended  by  having  the  accelerating  potential 
(3  kV)  turned  OFF  whenever  the  instrument  is  looking  into  the 
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Fig  . 


sun.  This  is  determined  by  a  small  photo-transistor  mounted 
behind  the  middle  slits  of  the  detectors. 

This  payload  has  only  one  operational  mode.  All  changes 
are  internal  and  only  ON  and  OFF  commands  are  required. 

Figure  2.6  is  the  electrical  block  diagram.  The  dc-dc  con¬ 
verter  supplies  power  through  the  doublers  to  the  plate  voltage 
stepper.  This  supplies  eight  voltage  steps  each  of  32-msec 
duration.  Each  voltage  step  is  monitored  to  the  telemetry. 
These  steps  are  supplied  to  each  set  of  analyzer  plates  as 
shown.  The  +3  kV  is  applied  to  the  channeltrons  of  each  ana¬ 
lyzer.  Collected  electrons  produce  a  cascade  effect  in  the 
channeltrons  and  the  resulting  pulse  is  amplified  and  then 
passed  to  a  logarithmic  accumulator  to  be  counted.  Each 
accumulator  produces  a  5-bit  mantissa  and  a  3-bit  exponent 
and  has  a  capacity  of  8032  counts.  Counts  are  accumulated 
for  the  32  msec  of  each  voltage  step,  at  the  end  of  which  the 
contents  of  all  four  accumulators  are  transferred  in  parallel 
to  the  output  shift  register.  The  accumulators  are  then 
zeroed,  the  deflection  step  voltages  are  incremented,  and  the 
accumulators  are  reenabled  for  counting  at  the  next  energy 
levels.  The  contents  of  the  shift  register  (32  bits)  are  read 
out  during  the  following  data  interval.  Since  there  are  eight 
voltage  steps  in  each  sensor  for  each  voltage  sweep,  a  sequence 
of  eight  32-msec  data  intervals  produces  two  complete  16-chan¬ 
nel  spectra  of  electrons  in  the  energy  ranges  of  50  cV  to  1  keV 
and  1  keV  to  20  keV  (one  spectra  for  each  of  the  two  look 
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directions).  The  data  words  from  each  analyzer  are  output  to 
Mainframe  91,  the  low  energy  analyzer  of  Detector  1. 

Figure  2.7  shows  the  CRL-251  enclosure.  Finally  Table 
2.2  gives  a  summary  of  the  instrument  characteristics. 

2.2.3  Summary  of  Operation 

The  P78-1  satellite  was  launched  on  February  24,  1979. 
"Quick  look"  data  indicate  that  the  electrostatic  analyzer  is 
operating  nominally.  Production  data  should  start  arriving 
at  AFGL  during  FY80. 
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Table  2.2 

ELECTROSTATIC  ANALYZER  (ESA) 
CRL-251 


Summary  of  Characteristics 


Particles  Detected 


Electrons 


Number  of  ESAs 


Four  (two  for  each  of  two-orthogonal 
look  directions) 


Detectors  Used 


2  Channeltron  Electron  Multipliers 
for  each  of  the  four  sets  of  ESA  plates 


Method  of  Energy  Analysis 


Voltage  Stepping  on  ESA  plates 


Number  of  Energy  Bins 


8  for  each  set  of  ESA  plates 


Energy  Bins  for  each  direction  Large  ESA  Small  ESA 

20.000kev  1.000  kev 

13.037  .652 

8.498  .425 

5.539  .277 

3.611  .181 

2.354  .113 

1.534  .077 

1.000  .050 


Geometric  Factor 


-5  *> 

Large  ESAs  :  4.85  X  10  ,cm  -ster 
Small  ESAs  :  1.30  X  10 


U  E/E 


Large  ESAs  :  4.2% 
Small  ESAs  :  6.7% 
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Data  Rato 


One  complete  spectrum  per  1/A  sec 
for  each  look  direction 


Dwell  Time  at  Each  Energy  Level  32  msec 


Digital  Data  Format  (8  bits  per  channel)X(8  energies  perESA)] 

X(2  ESAs  per  direction)X (2  directions)* 
=256  per  1/A  sec 


Note:  Each  8  bits  in  every  channel  are  in  logarithmic  form,  the  five  least  significant 
bits  being  the  mantissa  and  the  remaining  three  the  exponent.  This  number  is 
converted  to  decimal  form  according  to  the  following  relationship: 

N=  2y  (  x+32  )  -  33 

Where:  000  00000 <r'LSB 

y  x 


Weight 


A . 3A  lbs 


Sice 


6  in  X  6  in  X  A. 137  in 


Power  Dissipation 


.25  Watts 


+  28  Volts 
+  5  Volts 

ESA  plate  voltages 
Temperature  1 
Temperature  2 


Analog  Monitots 


3.  SIMULATION  OF  LEPS  RESPONSE  IN  ANY  GIVEN  PROTON 
BACKGROUND 

Considerable  amount  of  proton  data  were  obtained  by  the 
Low  Energy  Proton  Spectrometer  (LEPS)  instruments  flown  on 
board  the  S3-2  and  S3-3  satellites. 

However  not  all  of  the  data  were  useful  due  to  saturation 
and  energy  distortion  problems  which  were  encountered  in  areas 
where  large  proton  fluxes  exist.  An  attempt  was  made  at  simu¬ 
lating  the  LEPS  response  to  a  high  proton  flux  environment, 
with  mixed  success:  although  the  technique  used  was  accurate 
to  a  high  level  of  approximation,  it  could  not  be  made  to 
account  for  all  observed  artifacts  in  the  data. 

The  technique  used  in  the  simulation  of  the  LEPS  assumed 
that  the  proton  energy  spectra  as  observed  by  the  narrow 
angle  instrument  on  S3-3  was  sufficiently  accurate  to  do  a 
Monte  Carlo  type  partial  generation,  and  then  determine  the 
observed  spectrum  and  count  rates  in  the  wide  angle  instru¬ 
ment.  The  simulation  was  written  in  FORTRAN  language  on  the 
CDC6600  computer. 

3.1  Model  Description 

This  modeling  of  the  LEPS  instrument  is  divided  into 
two  parts.  The  first  part  characterizes  the  response  of  the 
instrument  in  a  given  angular  (solid  angle)  distribution  of 
a  proton  flux  for  a  given  orientation  of  the  instrument  with 
respect  to  this  distribution.  The  second  part  calculates 
the  response  of  the  instrument  to  the  above-mentioned 
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angular  distribution  for  any  given  energy  distributions. 

Two  programs  have  been  developed,  which  correspond  to  these 
two  parts.  The  first  program,  MODEL,  is  run  for  every 
angular  distribution  and  orientation  which  is  to  be  simu¬ 
lated.  For  each  run  of  program  MODEL,  various  energy  dis¬ 
tributions  can  be  then  run  by  means  of  the  second  program, 
ENERGY . 

As  the  programs  are  now  structured,  the  angular  distri¬ 
bution  used  is  a  so-called  pitch-angle  distribution  which  is 
of  the  form 


f(ft)  =  A^sinniJj 

where  A  is  the  normalization  constant  and  \p  is  the  angle 
between  the  particle  path  and  the  magnetic  field  vector  B. 
Program  ANGLE  asks  the  user  to  specify  both  the  'n'  in  the 
pitch  anfle  distribution  and  the  angle  between  B  and  the 

instrument  axis  z. 

The  energy  distribution  used  is  a  power  law  spectrum 
of  the  form 

f(0)  =  aee"n 

where  A  is  a  normalization  constant,  E  is  the  incident 
energy,  and  N  characterizes  the  power-law  spectrum. 

Other  analytic  functions  could  be  substituted  for  these 
functions  with  just  a  few  programming  changes.  Empirical 
data  (probably  in  the  energy  spectrum,  e.g.,  the  Vette 
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model)  could  also  be  substituted  with  a  few  more,  but  not 
substantial,  programming  changes.  Furthermore,  the  changes 
could  be  made  so  that  the  empirical  data  could  be  read  in. 

The  geometry  of  the  instrument  has  been  simplified  so 
as  to  be  more  easily  modeled  mathematically.  Because  the 
only  major  geometrical  difference  between  the  wide-  and 
narrow-angle  instruments  is  the  diameter  of  the  throat,  only 
this  change  need  be  made  to  the  dimensions  of  the  instrument 
used  in  the  program  in  order  to  have  the  program  model  one 
instrument  instead  of  the  other. 

Individual  ray  paths  are  selected  by  the  Monte  Carlo 
method  and  traced  geometrically  thru  the  instrument  to 
determine  the  path  length  thru  the  shielding,  the  path 
length  thru  the  first  detector,  and  whether  or  not  the 
second  detector  is  hit  (i.e.,  whether  or  not  the  anti- 
coincidence  circuit  is  triggered).  The  path  length  thru 
various  materials  used  in  the  shielding  has  been  reduced  to 
an  effective  path  length  in  silicon.  By  means  of  the  appro¬ 
priate  angular  distribution  and  orientation  of  the  instru¬ 
ment  axis,  program  MODEL  computes  and  stores  this  informa¬ 
tion,  which  is  then  put  on  a  disk  memory. 

Program  ENERGY  takes  the  information  supplied  by  program 
MODEL  and  calculates  the  response  of  the  instrument  for 
various  energy  spectra. 

The  assumptions  inherent  in  this  model  are  the  following 

1)  The  protons  move  in  straight  lines  as  they  pass  thru 


the  shielding.  (Note  that  any  other  sub-atomic  particle  could 
be  modeled  as  long  as  (a)  it  travels  in  a  straight  line  thru 
matter,  and  (b)  it  has  a  known  range-energy  distribution  with¬ 
out  appreciable  statistical  fluctuation.) 

2)  There  are  no  saturation  effects  taken  into  account. 

3)  Any  effects  where  the  back  detector  only  is  hit  are 
not  taken  into  account. 

4)  Only  particles  from  the  'forward  hemisphere'  of  the 
instrument  are  taken  into  account.  (Note  that  the  'back 
hemisphere’  could  be  added,  but  this  would  imply  a  knowledge 
of  the  effective  shielding  of  the  rest  of  the  satellite.) 

5)  The  spectrum  f(E,  )  is  assumed  to  be  separable  as 

f(E,fl)  =  f  t (E)  f2(n) 

i.e.,  the  particular  angular  dependence  (or  non - i sotropy)  is 
independent  of  the  energy  of  the  particle. 

3.2  Modelling  the  Instrument  Geometry 

To  model  the  LEPS  on  the  computer  it  is  advisable  to 
reduce  the  instrument's  relatively  complex  geometry  to  a  sim¬ 
plified  form  which  is  more  tractable  to  being  modeled  in 
FORTRAN  on  the  computer.  A  picture  of  the  LEPS  is  shown  in 
Figure  1.  With  the  exception  of  the  magnets  placed  on  either 
side  of  the  throat,  the  instrument  is  cyl indrical ly  symmetri¬ 
cal.  A  cross-section  of  the  LEPS  is  shown  in  Figure  2.1.  In 
this  cross-sectional  view  the  geometry  of  the  instrument  has 
been  simplified,  and  it  is  quite  close  to  what  is  actually 


modeled  on  the  computer. 

To  calculate  the  path  length  thri  the  shielding  of  a 
particular  ray,  it  is  necessary  to  know  what  surfaces  the 
ray  has  penetrated  and  what  material(s)  were  penetrated.  A1 
particles  are  assumed  to  hit  the  front  detector.  (If  a 
particle  does  not  hit  the  front  detector,  it  is  not  counted 
anyway.)  The  rays  are  therefore  grouped  into  categories 
depending  on  which  surfaces  it  has  penetrated,  and  the  ray 
path  is  then  calculated  from  the  ray-path  parameters  based 
on  the  category  into  which  it  fits. 

The  first  separation  is  based  on  whether  the  front 
detector  is  hit  at  a  point  inside  or  outside  the  throat. 
Within  each  of  these  categories  the  rays  are  separated  based 
on  which  surfaces  have  been  penetrated. 

In  the  computer  program  the  surfaces  which  are  pene¬ 
trated  are  determined  from  the  parameters  of  the  ray  path. 
Thru  each  part  of  the  shielding  the  path  length  is  then  cal¬ 
culated  in  cm  based  on  these  ray-path  parameters. 

3.3  Range-Energy  Considerations 

Since  the  instrument  shielding  is  composed  of  different 
materials,  it  is  convenient  to  reduce  the  path  length  thru  a 
particular  material  to  an  equivalent  path  length  in  any  one 
given  material.  Jilicon  has  been  chosen  since  this  is  what 
the  detector  consists  of.  Hence  the  path  length  in  cm  thru 
each  piece  of  the  shielding  is  converted  to  an  equivalent 
path  length  in  silicon  for  that  piece.  The  path  lengths  are 
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then  added  to  produce  a  total  equivalent  path  length  thru  the 
shielding  in  equivalent  cm  of  silicon  for  that  particular  ray 
path. 

The  reason  that  a  path  length  thru  a  certain  material 
can  be  reduced  to  an  equivalent  path  length  in  silicon  is 

based  on  the  empirical  fact  that  when  the  range,  expressed  in 

2 

gm/cm  ,  is  plotted  versus  the  energy  of  the  incident  particle, 
for  protons  in  MeV,  the  curves  for  all  the  materials  in 
question  are  the  same  within  a  relatively  small  error.  The 
path  length  thru  each  material  in  cm  is  multiplied  by  an 
appropriate  conversion  factor,  this  factor  being  (almost)  in¬ 
dependent  of  energy. 

3.4  Ray-Paths  —  Selection  Process  by  Monte  Carlo  Method 

All  particles  are  assumed  to  pass  thru  the  front  detector 
since,  if  a  particle  does  not  hit  the  front  detector,  it  is 
not  counted.  Particles  entering  the  back  surface  of  the  front 
detector  are  discounted.  In  reality  the  effective  shielding 
of  the  satellite  would  have  to  be  taken  into  consideration, 
and  this  is  not  at  any  rate  known  to  any  degree  of  accuracy. 
Particles  entering  the  "sides"  of  the  detector,  producing 
what  are  called  'side  tracks',  have  been  considered  but  are 
not  dealt  with  in  this  report.  Hence  only  particles  which 
pass  thru  tho  front  surface  of  the  front  detector  are  con- 
s i dered . 

Hence,  a  ray  path  is  chosen  by  selection,  a  point  on 
the  surface  of  the  front  detector  and  a  direction  in  space  by 


the  Monte  Carlo  method.  The  point  on  the  surface  is  chosen 
at  random  based  on  the  assumption  that  the  flux  is  uniform 
across  any  given  surface.  The  direction  is  chosen  at  random 
based  on  an  isotropic  flux,  i.e.,  the  flux  is  proportional 
to  the  solid  angle  as  seen  by  the  point  on  the  detector  sur¬ 
face  which  is  hit.  The  effective  cross-section  of  the  detec¬ 
tor  surface  and  any  deviations  from  an  isotropic  distribution 
are  put  in  as  weighting  factors  as  explained  below. 

3.5  Coordinate  Systems 

Let  the  origin  of  any  coordinate  system  be  at  the  center 
of  the  first  detector  with  the  z-axis  the  plane  of  the  detec¬ 
tor.  Because  of  the  cylindrical  symmetry  of  the  model,  we 
need  not  specify  the  orientation  of  the  x  and  y  axes.  (How¬ 
ever,  if  for  instance  the  magnets  are  added  to  the  model, 
these  axes  would  have  to  be  defined  relative  to  the  posi¬ 
tions  of  the  magnets.) 

We  can  use  coordinate  systems  wh  ch  are: 

Cartesian  x,y,z 

Cy 1 indrical  r  ,  0  ,  z 

Spherical  p , a , B 

These  conventions  are  used  throughout. 

3.6  Random  Number  Selection 

To  select  a  parameter  x  randomly,  where  xt  4  x  ^  x2, 


weighted  by  a  function  f(x)  we  have 


dr'  =  f(x)dx 


r  x 

f  dr'  =  A  /  f(x')dx' 


(i: 


where  r  is  the  random  number  selected,  0  <  r  <  1,  and  A  is  a 
normalization  factor  defined  by 


1  =  A  /  f (x ' ) dx  ' 

xi 


and 


A 


1 


f (x ' ) dx  1 


C  imbining  (1)  and  (2)  we  have, 


(2; 


r 


x 

/  f(x')dx' 

xi 


I 

X  ■ 


f  (x  '  ) dx  ' 


(3) 


3.7  Selection  of  Ray  Path  Pa r ameters  by  Means  of  the 
Monte  Carlo  Method 

The  ray  paths  are  chosen  by  random  numbers.  Two  random 
cumbers  can  be  used  to  specify  a  point  on  the  front  detector 
1 r | , 0 i )  in  cylindrical  coordinates,  and  two  more  to  specify 
;  direction  in  space  (uj.gj)  in  spherical  coordinates.  Let 
the  random  numbers  be  j  where  i  is  the  number  of  the  ray 
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path  and  j  is  1-4,  the  particular  random  number  of  the  set  1 
We  then  have 

1)  for  r 


dft . ,  =  r .  dr . 

1 1  l  l 


f  r .  dr  2 

s  =  o _ =  !i 

11  r0  rD 
/  ri  dr 

0 


where  rp  =  radius  of  front  detector 


ri  =  rd/fiil 


2)  for  6 


= 


/  d6, 

f  d  0  i 
0 


0j  =  2irni2 


3)  likewise  for  aj  we  have 


04  =  2tt«  i  3 


4)  for 


cos0  dgj 


/  cosg.  dg. 

-  TT  /  2  1  1  sing.  +  1 

°i4  ~  Y/l  "  “ 

/  cosg.  dg.  1 

-•nr/ 2  1  1 


sing.  =  2Q. .  .  -  1 
l  1 4 


3.8  Ray-Path  Equation 


We  now  have  a  point  on  the  surface  of  the  detector 


which  we  denote  by  and  a  direction  in  space  d. 
Hence  the  equation  of  the  ray  path  is 


P  =  P0  +  sd 

where  p  is  the  position  of  a  given  point  with  respect  to 
the  origin 

Pq  is  the  position  of  an  initial  point  on  the  detector 
with  respect  to  the  origin 

d  is  a  unit  vector  defining  the  direction  in  space 
s  is  the  path  length  from  the  detector  surface 
in  terms  of  r,  0,  a,  and  we  have  (note  that  the  subscript 
i  has  been  dropped)  : 

A  A  A 

p0  =  x  r0cose0  +  y  rQsin0o  +  zzQ 

A  A  A 

d  =  x  cosacosg  +  y  sinacosg  +  z  sing 
and  therefore 


p  =  x(rQ  cos0q  +  s  cosacosg) 

A 

+  y (rQ  sin0Q  +  s  sinacosg) 
+  z(s  sing  +  z q ) 


the  path  length  can  be  expressed  in  terms  of  x,  y, 
or  z  as 

x  -  ro  cos60  _  y  -  r0  sin90  _  z  -  20 

cosacosB  sinacosB  sing 

where  the  subscript  'O'  is  used  to  note  the  value  of  a 
parameter  on  the  detector  surface.  The  equations  for 
x,  y,  and  z  are: 

x  =  Tq  co s 0 q  +  s  cosacosB 
y  =  r„  sin9  +  s  sinacosB 

u  o 

z  =  zQ  +  s  sinB 

the  path  length  (s2~Sj)  between  two  points  1  and  2 
is  therefore: 


cosacosB  sinacosB  sin0 
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where  A  =  cos  8 


B  =  2rQ  cos8cos (0Q-a) 

2  2 

c  =  Tn  -  r 
0  s 


The  intersection  (if  there  is  one)  of  a  ray  with 
a  cone  can  be  found  in  similar  fashion. 


3.9  Calculation  of  the  Effective  Path  Length  Thru  the  Shielding 

3.9.1  Intersection  of  the  Ray  Path  with  a  Given  Surface 

To  solve  for  the  intersection  of  the  ray  path  with  a  given 
surface,  we  take  the  equation  of  the  surface  and  plug  it  into 
the  equation  for  p.  For  this  instrument  the  geometrical 
figures  of  interest  are  the  cylinder,  plane  surface  at  a  given 
z,  and  the  cone.  This  is  part  of  the  rationale  for  the 
simplification  of  the  instrument  geometry  explained  above. 

The  equation  for  p  is  then  solved  to  determine  s,  the 
distance  from  the  surface  of  the  detector  at  which  the  ray 
path  intersects  the  surface  of  interest.  An  example  is  here 
given  for  a  surface  z  =  z$,  i.e.,  one  which  is  parallel  to 
the  surface  of  the  detector,  and  r  =  r  ,  i.e.,  a  cylindrical 
surface  which  is  concentric  with  the  detector. 

3.9.2  Calculation  of  the  Effective  Path  Length  Thru 
the  Shielding 

Knowing  the  points  1  and  2  at  which  the  ray  intersects 
the  two  appropriate  surfaces  which  delimit  a  particular  piece 
of  the  shielding,  one  can  then  calculate  the  path  length  in 
cm  through  this  piece  by 

As  =  s „  -  s , 
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This  is  then  converted  to  equivalent  cm  in  silicon  s’  by 
multiplication  by  the  appropriate  factor  f 


As*  =  fAs 

The  sum  of  the  effective  path  lengths  Rs  is  then 

R  =  E  As '  =  E  f .  s.' 

s  .11 

1 

If  the  particle  comes  straight  down  the  throat,  no 
shielding  is  penetrated,  Rg  =  0,  and  the  ray  is  termed 
a  'legitimate'  ray.  Conversely,  for  Rg  >  0,  the  ray  is 
termed  an  'illegitimate'  ray. 


3.10  Calculation  of  the  Path  Length  Thru  the  Front  Dei e c t o r 
The  path  length  thru  the  front  detector  R^  is  then 
calculated  by  the  formula 


R 


D 


sing 


where  t^  is  the  thickness  of  the  front  detector.  Note  that 
this  calculation  assumes  that  the  ray  passes  out  at  the  front 
detector  thru  its  back  surface  and  not  thru  the  side.  This 
is  a  reasonable  approximation  to  the  actual  case  since  not 
many  rays  will  pass  thru  the  side.  This  effect  is,  however, 
taken  into  account  when  the  sides  of  the  detector  are  considered 
and  is  not  treated  in  this  report. 

3.11  Weighting  Factor  and  Rg  vs  RR  Matrices 

These  two  path  lengths  and  whether  the  anticoincidence 
is  triggered  provide  the  necessary  information  to  determine 
the  response  of  the  instrument  to  a  particle  with  the  given 
particle  path  with  a  given  energy.  However,  program  MODEL 


takes  these  quantities  and  enters  a  number  in  one  of  two 
matrices,  one  for  which  the  anticoincidence  is  triggered 
and  one  for  which  is  not.  The  element  of  the  appropriate 
matrix  to  which  the  number  is  added  is  determined  by  the 
path  length  thru  the  shielding  Rs>  which  determines  the  row 
of  the  matrix,  and  the  path  length  thru  the  first  detector 
Rp,  which  determines  the  column  of  the  matrix.  The  number 
which  is  actually  entered  in  this  matrix  element  is 
determined  by  the  angular  distribution  itself  and  by  the 
effective  cross  section  of  the  front  detector. 

The  upper  and  lower  limits  for  Rg  and  Rp  are  predetermined, 
and  the  'boxes'  of  the  matrices  are  evenly  spaced  logarith¬ 
mically.  For  example,  if  R&  and  Rp  are  a  pair  of  path  lengths 
we  have 


log  R, 


log  (R  )  -log  (r  ) 

B  s'max  BV  s  mm 


(7  -1) 

max 


log  R 


<  I  < 


log(R)  -log(R)  . 

6  s 'max  s 'min 


max 


log  R 


D 


log ( R  )  -log(Rn) 

6  D  max  BV  D'min 


(J  -1)  <  J  < 
max  '  — 


log  R 


D 


log(Rn)  -logfRpJ  . 

6  D'max  6  D  min 


J 

max 


where  I  and  J  denote  the  matrix  element,  and  I  and  J 

max  max 

are  the  number  of  rows  and  columns  in  the  matrix. 

The  selected  ray  is  chosen  with  no  regard  as  to  the 

effective  cross  section  of  the  detector  surface  or  the  angular 

distribution  in  question.  These  effects  are  added  by  means  of 

a  weighting  factor.  The  effective  area  of  the  detector  is  cos3, 

and  the  angular  distribution  function  f  (SI)  provides  its  weighting 

2 

factor.  In  the  case  where  f  (ft)  =  sin  ip  we  have  for  a  weighting 
2 

factor  cosgsin  tp  where  ip  =  \p(a,B)  is  known.  This  function  ip 
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will  depend  on  ipQ ,  the  angle  between  the  inst  rument  axis  z 
and  the  magnetic  field  B. 

This  weighting  factor  is  then  added  to  the  appropriate 
element  I,  J,  calculated  above,  of  one  of  the  two  matrices, 
depending  on  whether  or  not  the  ant i - co inc idence  circuit  is 
triggered.  The  program  then  goes  back  to  select  the  next 
ray.  The  two  matrices  are  thus  built  up  as  the  requisite 
number  of  rays  is  processed.  (10,000  rays  seems  to  give 
good  results,  and  this  takes  5  seconds  of  computer  time. 

For  better  accuracy  100,000  rays  are  recommended,  which 
would  take  50  seconds.) 

If  the  ray  is  a  legitimate  ray  the  weighting  factor 
is  added  to  the  legitimate  ray  count. 

Note  that  in  Program  MODEL  the  ant i - coinc idenc e  matrix 
is  used  as  long  as  the  second  detector  is  hit.  There  are 
instances  whrere  the  anticoincidence  circuit  is  not  triggered 
even  though  the  particle  hits  the  second  detector.  However, 
effects  such  as  these  are  taken  into  account  in  Program 
ENERGY. 

This  concludes  the  calculation  done  by  the  first  Program 
MODEL .  The  matrices  and  other  information  are  printed  on  the 
line  printer  and  are  stored  on  a  permanent  file  of  the  disc 
to  be  processed  by  the  second  Program,  ENERGY .  Given  the 
particular  instrument  to  be  modeled,  wide  or  narrow - ang 1 e , 
and  given  the  pitih-angle  distribution  and  the  angle  of  the 
instrument  axis  with  the  magnetic  field,  MODEL  provides  the 
number  of  legitimate  particles  and  histograms  of  path  length 
thru  the  shielding  versus  path  length  thru  the  detector  for 
illegitimate  particles  for  both  the  coincidence  and  anti- 


coincidence  mode. 


3.12  Construction  of  the  Joint  Distribution  of  R  and  R^ 
_ s _ D 

The  input  to  program  ENERGY  is  the  number  of  legitimate 
particles  and  two  matrices  which  are  number  densities  as  a 
function  of  Rg  and  RQ  for  the  anticoincidence  (AC)  mode  and 
the  non-anticoincidence  (NAC)  modes,  respectively.  If 
and  represent  these  two  matrices,  respectively,  and 

L  represents  the  number  of  legitimate  particles,  the  total 
number  of  particles  is 

N  =  L  +  E  N.  .  .  +  E  N_ .  . 

P  ^  •  I13  •  •  2ii 

1.3  1,3  J 

Dividing  the  elements  of  each  of  these  matrices  by  N  , 
we  have  an  approximation  to  a  joint  probability  distribution. 
Defining  the  mode  m  such  that 

m  =  0  Legitimate  particle 

m  =  1  Illegitimate  particle  with  anticoincidence 
m  =  2  Illegitimate  particle  with  no  anticoincidence 

we  can  express  this  distribution  as 

P (0 , 0 , 0)  =  Jp 


Rs(i+1)  RD(i+l) 


p(Rs,RD,m)  dRsdRD 


N  .  . 
mi  j 


m  =  1,2 


Rs(i)  R„(j) 


or 


P(Rs,RD’n)  AR  AR 


N  . 
mi  j 


D  "p 


where  AR^AR^is  the  'area'  of  a  matrix  element.  Note  that 
these  'areas'  will  vary  from  element  to  element  since  they 
are  set  up  logarithmically. 


*  •* 
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3.13  Determination  of  the  Energy  Dropped  in  the  First  Detector 

The  joint  probability  distribution  p(Rs,RD,m)  is  a  function 
of  Rg  and  R^.  It  is  necessary  to  find  the  energy  dropped  in 
the  first  detector  for  a  given  incident  energy  Eq  as  a  function 
of  Rs  and  Rp.  Then  given  Eq  we  can  determine  the  conditional 
distribution  pCE^/Eq),  i.e.,  the  probability  of  an  energy  E^ 
being  deposited  in  the  front  detector,  given  an  incident 
energy  EQ . 

This  is  done  through  application  of  the  range-energy 
relation  (derived  empirically  from  Barkas  and  Berger) 


where  R  is  the  range  in  silicon  for  a  proton  of  incident 
energy  Eq,  and  B  and  M  are  constants  derived  from  the  fit 
to  the  empirical  data.  Since  Rg  is  understood  to  be 
equivalent  cm  of  silicon,  and  R^  is  the  range  through  a 
silicon  detector,  this  formula  is  all  that  is  needed. 

By  algebraic  manipulation  we  find  that  the  energy  E^ 
of  the  particle  after  it  leaves  the  shielding  is 


The  energy  the  particle  has  lost  in  the  shielding  is 


E  ' 


1/M 


It  is  next  determined  whether  or  not  the  particle  is  stopped 


in  the  first  detector.  If  the  relation 


is  true,  the  particle  is  stopped  in  the  first  detector  and 
we  have 


=  E '  =  E , 


For  Rq  <  BE^M,  we  must  find  Eq  in  a  manner  similar  to  the  way 
E'  is  found,  and  we  have 


E  p.  =  E 
D  x 


,M 


1/M 


Since  the  calculation  of  Eq  is  not  straightforward,  and 
it  is  not  desirable  to  use  the  analytic  expression  for  Eq, 
it  becomes  necessary  to  choose  a  number  of  values  for  Eq 
which  adequately  cover  the  energies  of  interest  at  the  in¬ 
coming  flux.  Then  for  each  value  of  Eq  a  table  is  made  of 

ed  =  ed^rs,rd^e=ed 


where  the  values  of  Rs  and  Rq  in  the  table  correspond  to 
those  input  from  program  ANGLE . 

At  present  15  values  of  Eq  are  calculated.  These  values 
span  the  interval  from  10  to  500  MeV  in  a  quasi  - 1 ogarithmic 
manner.  The  actual  values  are  10,  20,  30,  40,  50,  60,  70, 

80,  90,  100,  125,  150,  175,  200,  300,  400,  and  500  MeV. 

3.14  Calculation  of  the  Conditional  Distribution 

For  each  value  of  Eq  chosen,  the  table  of  E  q  ( R  s  ,  R  q  )  is 
scanned,  and  for  all  entries  of  Eq,  which  fall  within  the 
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is  added 


channel  width  of  interest,  the  value  of  N^j  or  N2ij 
to  make  up  the  conditional  distribution  as  explained  below. 

We  here  define  ’n’  as  the  channel  number  and  en  as  the  thresh¬ 
old  energy  of  the  nth  channel.  (The  cutoff  energy  of  the  top, 
or  12th  channel  is  6  MeV) .  We  are  then  calculating  a  function 
Fc(EQ,n),  which  corresponds  to  a  theoretical  conditional  prob¬ 
ability  distribution  of  Eq  given  Eq,  which  we  will  write 
4>(Eq/Eq).  We  then  have 


Fc(E0,n)  =  Z  N2  i  j  +  Z  Nj  ifj 
i , j  i,  j 

e  ED  En+1  ED  =  Eex 

The  first  term  in  this  equation  counts  rays  for  which  the 
anticoincidence  circuit  was  not  triggered  due  to  the  fact  that 
the  ray  misses  the  second  detector:  this  dictates  the  use  of 
the  matrix  N2jj.  The  second  term  counts  rays  for  which  the 
anticoincidence  circuit  was  not  triggered  due  to  the  fact  that 
the  ray  was  stopped  in  the  front  detector  -  hence  the  use  of 

Nlij- 

Relating  Fc(EQ,n)  to  the  theoretical  distribution 
(J>  (Eq/Eq)  we  have 


F 


c 


E 


0  > 


N 


n) 


en+  1 

/ 


d>(ED/Eo)  dE0 


We  can  calculate  an  approximation  to  tp  (E^/Eq)  by 


<M|:.|) /lip) 


Fc(E0>n) 

N|>  n+  1  ‘  1  n  1 


Relating  it  to  the  joint  probability  distribution  we  have 


^c(^0,n)  En+1 

~  =  f  dE  ff  p(R  R  2)  dR  dR 

en  e  .  s  u 

n  n  0  n+ 1 


Ni 


'n+  1 


+  /  dEn  ff  p(Rs,RD,l)  dRs  dRD 


En  =  E 
0  ex 


3.15  The  Incident  Spectrum 

Now  the  conditional  distribution  <p  (  Eg/Eg)  is  related  by 
definition  to  the  joint  distribution  ffE^.E^)  by  the  equa¬ 
tion 

f(ED,E0)  -  <KEd/E0)  f0(E0) 

where  fgCEg)  is  the  marginal  distribution  of  Eg.  Physically 
fg(E0)  is  equivalent  to  the  incident  spectrum,  i.e.,  the 
distribution  of  incident  energies  of  the  protons,  as  long  as 
the  following  normalization  is  maintained. 

f  fg (Eg)  dEg  =  1 

where  Ej  and  E ^  arc  the  bounds  of  the  energy  range  of  interest 
E.  is  set  to  Ei,  the  threshold  Ei  of  the  first  channel 


- - 


(  =  100  keV)  ,  and  E2  is  set  to  the  maximum  energy  for  which  there 
are  still  appreciable  counts  in  the  incident  spectrum.  We  have 
chosen  E2  =  500  MeV . 

In  the  program  as  it  now  stands,  a  power  law  distribution 
is  used  of  the  form 

fo(Eo>  '  sboN 

where  N  characterizes  the  distribution  and  A  is  a  constant 
found  by  solving  the  normalization  condition  explained  above. 
Power  law  distributions  with  N=l,2,4,8  are  used  in  the  program. 
By  means  of  fairly  simple  programming  changes,  this  function 
could  be  changed  to  another  empirical  function,  or  empirical 
data  with  a  simple  interpolation  scheme  could  also  be  substi¬ 
tuted. 

The  number  of  legitimate  particles  Cj(n)  in  a  particular 
channel  can  now  be  calculated  by 

ci<">  ■  i  ^n*'  fo(E>  dE 

1  e 

n 

3.16  The  Background  Count 

The  background  which  will  be  counted  by  the  LEPS  is  equal 
to  the  marginal  distribution  of  the  energy  E  p  deposited  in  the 
first  detector,  fp(Ep).  This  can  be  found  by  integrating  the 
joint  distribution  of  Ep  and  Hq  over  the  whole  energy  spectrum, 

W  ■  {E"aX  f<ED’E0>dE0 
E  . 
mm 
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To  get  the  background  counts  in  a  given  channel  Cg(n),  we  must 
integrate  with  respect  to  E^  over  that  channel's  range  of 


energy,  and  we  have 


CB<">  '  !  f0(E0)dE0 

e„ 


=  /  n+1dED  /  maX  f(ED,EQ)dE0 

e  E 

n  min 


In  the  program  we  calculate  an  approximation  to  Cg(n)  by 

CR(n)  =  i  E  F  (E . ,  n)  fn(E.)  (E.  .  -  E.) 

B  N  _  c  v  1  ’  '  0  1  l+l  r 

Ei 

where  E^  are  the  incident  proton  energies. 

For  each  channel,  the  total  count  C^(n)  is  the  sum  of 
the  background  and  the  legitimate  particles 

CT(n)  =  Cg  (n)  +  CL(n) 

the  fraction  p(n)  of  legitimate  particles  in  the  total  count 
is  then 

Cl>) 

P(n)  =  CB(n)  +  C  L (n ) 

when  looking  at  real  data,  it  is  this  quantity  to  which  the 
data  would  be  compared. 


3.17  Conclusions 

The  modeling  which  has  been  developed  will  be  used  to 
interpret  data  which  have  been  taken  by  existing  I, EPS  instru¬ 
ments  in  orbit.  Spectra  measured  by  the  LEPS  will  be  input  to 


the  programs,  and  the  programs  will  produce  the  output  which 
the  respective  spectra  should  have  generated.  Comparison  will 
reveal  whether  the  data  actually  taken  were  good. 

Another  task  that  is  proposed  is  to  input  data  from  an 
existing  model  of  the  magnetosphere,  such  as  the  Vette  model. 
The  programs  will  provide  output  telling  what  the  instrument 
should  measure  in  the  given  environment.  This  output  will 
then  be  compared  with  actual  data  taken  by  the  LEPS.  In  this 
way  the  modeling  of  LEPS  in  conjunction  with  the  LEPS  data 
can  be  used  to  test  existing  models  of  the  magnetosphere. 

Further  modifications  can  also  be  made  to  take  into 
account  electronic  characteristics  of  the  instruments.  For 
instance,  the  program  could  be  programmed  to  take  into 
account  any  dead-time  effects.  This  type  of  extension  to  the 
modeling  would  possibly  require  some  laboratory  work  to 
exactly  define  the  response  of  the  instrument  to  a  particlar 
effect.  Once  this  effect  has  been  sufficiently  defined,  it 
would  then  be  modeled  mathematically  and  coded  into  the 


programs . 
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4.  DATA  ANALYSIS 

4.1  Data  Handling  Software  Development 

An  interest  was  expressed  in  scanning  the  data  bases  of 
the  S3-2  and  S3-3  satellites  for  evidence  of  solar  activity. 
In  order  to  accomplish  this,  some  restructuring  of  the  data 
files  was  necessary  to  allow  for  analysis  of  a  full  360° 
pitch  angle  distribution  and  remove  an  artificially  low 
upper-limit  on  the  "L-shell"  calculations.  This  restruc¬ 
turing  was  done  and  one  tape  containing  25  orbits  of  data 
from  the  S3-2  satellite  was  reprocessed  to  serve  as  input 
for  development  of  the  data  analysis  software. 

One  initial  goal  of  the  S3-2  and  S3-3  data  analysis 
projects  was  to  generate  a  complete  graphical  documentation 
of  the  data  bases  through  calcomp-type  CRT  plots.  This  goal 
was  abandoned  when  it  became  obvious  that  the  existing  CRT 
facilities  at  AFGL  were  inadequate.  However,  the  CRT/ 
Microfiche  equipment  recently  installed  at  AFGL  does  have 
the  capability  we  need  for  the  project.  Most  of  the  re¬ 
quired  software  was  developed  before  the  original  abandon¬ 
ment  of  the  project.  However,  our  plotting  software  had 
to  be  updated  to  provide  options  for  the  additional  capabil¬ 
ities  of  the  new  CRT  facility.  Also  the  original  decision 
to  fold  the  pitch  angle  dependence  of  the  data  base  into  a 
180°  system  proved  to  be  a  mistake  since  it  did  not  allow 
for  adequate  visualization  of  the  polar  region  activity. 

Specifically,  updates  were  required  in  the  following' 
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data  reduction/data  analysis  programs. 

a)  S2BLAF/S3BLAC .  These  programs  add  the  magnetic  field 
and  pitch  angle  coordinates  to  the  data  files.  It  was  neces¬ 
sary  to  begin  reprocessing  the  data  bases  at  this  stage  be¬ 
cause  the  previous  cutoff  at  "L"  =  10.00  proved  to  be  too  low 
for  analysis  of  the  polar  region  data. 

b)  S2SBFK/S3SACD .  These  programs  separate  the  data  into 
sequential  groups  relative  to  predetermined  "L"  and 
boundaries.  Alterations  were  necessary  to  incorporate  the 
new  360°  data  formats. 

c)  S2PLAA/S3PLAA .  These  programs  process  the  output 
from  the  above  stage  and  generate  an  intermediate  file  for 
input  to  our  plotting  software.  An  update  was  necessary  in 
order  to  incorporate  the  360°  data  format  and  to  maintain 
compatabi 1 i ty  with  our  plotting  software. 

d)  DCRVBH.  This  program  and  its  associated  subroutines 
is  our  basic  plot  generating  software.  Several  revisions 
were  made  to  it  relative  to  the  new  capabilities  provided  by 
the  recently  installed  CRT/Microfiche  facility  at  AFGL.  Also 
a  major  revision  was  made  to  enhance  its  frame  blocking  cap¬ 
abilities  so  that  it  is  now  possible  to  combine  up  to  15 
plots  into  one  frame  in  either  a  1  by  <15,  2  by  <5,  or  3  by 
<5  format.  The  previous  version  had  only  the  1  by  <15  cap¬ 
ability. 

The  above  revisions  were  complete  for  the  S3-2  satellite 
and  the  test  tape  was  run  giving  satisfactory  results. 
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However,  the  decision  was  made  to  proceed  directly  to  the  S3-3 
data  set,  incorporating  the  ideas  developed  for  S3-2  into  the 
S3-3  graphics  package. 

4.2  S3-2  Proton  Data  Analysis 

4.2.1  Goal 

The  Low  Energy  Proton  Spectrometer  (LEPS)  on  board  the 
S3-2  polar  orbiting  satellite  provides  statistically  signifi¬ 
cant  proton  data  in  the  100-200  keV  energy  range  over  the 
polar  cap  regions.  These  data  are  available  over  a  six  month 
interval  from  December  1975  to  May  1976.  The  angular  resolu¬ 
tion  of  the  pitch  angle  distribution  is  fifteen  degrees,  and 
thus  provides  valuable  details  in  these  observations.  The 
purpose  of  this  investigation  was  to  study  the  different 
proton  pitch  angle  distributions  and  their  relation  to  the 
solar  activity,  magnetic  activity,  and  possibly  to  the  geo¬ 
graphic  longitude  in  connection  with  the  South  .Atlantic 
anomaly . 

4.2.2  Method 

Initially,  a  graphical  display  of  the  raw  data  for  the 
100  or  200  keV  protons  offered  a  "quick  look"  at  the  data  for 
particular  polar  passes.  Figures  4.1  and  4.2  show  plots  of 
the  count  rates  vs.  time  (in  seconds)  for  typical  passes. 
Varying  flux  distribution  patterns  are  seen  as  the  satellite 
moves  from  the  trapped  regions  where  the  flux  is  modulated 
by  the  sixteen  second  spin  period  of  the  satellite 
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(cf.  Fig.  4.3),  to  a  region  of  quas i - trapped  particles  which 
exhibit  a  flat  distribution  outside  the  loss  cone  (cf.  Fig. 
4.4).  Over  the  polar  regions  (A  >  70°),  we  observe  a  peak, 
at  approximately  | 120° |  in  the  pitch  angle  distribution 
(cf.  Fig.  4.6  to  4.9).  Figure  4.5  shows  a  typical  polar 
pass  in  invariant  latitude,  magnetic  local  time  coordinates. 

Plots  exemplified  by  Figures  4.1  and  4.2  enable  us  to 
set  the  spatial  limits  over  which  data  can  be  superimposed 
and  averaged  (quasi-trapped ,  polar  magnetic  night  or  day). 
These  averagings,  over  a  given  AL,  were  made  over  time 
intervals  of  approximately  one  hundred  seconds. 

For  the  purpose  of  uncovering  any  correlation  of  proton 
flux  with  solar  activity,  the  active  period  of  March  26  to 
April  1,  1976  was  selected. 

4.2.3  Results 

The  results  of  this  investigation  were  presented  at  the 
1978  Fall  AGU  Meeting  (cf.  EOS,  59^,  1  17  2  ,  1978).  During  the 
above-mentioned  period  of  high  solar  activity,  over  the  polar 
region  above  A  =  76°,  a  peak  repeatedly  appeared  in  the  pitch 
angle  distribution  at  approximately  -120°  (cf.  Figs.  4.6  to 
4.9).  This  peak  was  not  present  during  times  of  quiet  sun 
and  quiet  magnetic  field  on  April  17,  1976,  as  shown  in 

Figures  4.10  and  4  11.  On  this  date  the  count  rates  were 
very  low  (M  sec  ^)  for  all  pitch  angles.  This  fact  lead 
us  to  believe  that  there  was  no  solar  UV  contamination  of 
the  observations.  Also,  because  of  the  efficiency  of  the 
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sweeping  magnets  in  the  LEPS,  we  felt  confident  that  there  was  virtually 
no  electron  contamination  in  the  proton  data.  Further  analysis  of  the 
data  for  other  times  in  the  six  month  period  is  projected  in  order  to 
ascertain  the  possible  correlation  which  seems  to  exist  between  the 
peculiar  pitch  angle  distributions,  as  observed  from  March  26  -  April  1, 
1976,  and  the  solar  and  magnetic  activities,  and  possibly  with  the 
geographic  longitude. 

Documentation,  overhauling,  debugging,  and  testing  are  proceeding 
on  the  LEPS  Monte  Carlo  simulation  program.  Major  changes  are  being 
made  which  will:  1)  make  the  program  easy  to  use  for  those  not  familiar 
with  computers  and  2)  make  the  simulation  applicable  to  instruments 
other  than  the  LEPS. 

Whereas  before,  the  specific  geometry  of  the  LEPS  was  actually 
written  into  the  Fortran  code,  the  program  will  now  be  independent  of 
the  particular  dimensions  and  geometry  of  the  instrument,  as  long  as 
the  geometry  of  the  instrument  is  'similar'  to  that  which  the  program 
can  accept.  The  program  can  then  be  used  on  instruments  of  somewhat 
different  geometry  without  changing  any  of  the  Fortran  code. 

Debugging  and  testing  of  the  program  are  proceeding  concurrently 
with  the  above  work. 
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-4— f-f-U+H 


i  NORTH ) 


PITCH  ANGLE 

-2  LAPE  HULA  1589  3/27/76  112  SECONDS  DATA 

11/17/78  PLOT  NO.  9  UT=19.23  LT=11.35  MT=9.15 

®  *♦  -  5  MEV  &  1.2  MEV  +  .2  MEV  X  .1  MEV  $  75.  MEV 

0  37.5  MEV  X  19  MEV  Z  10.5  MEV  Y  6.8  MEV 


(  NORTH  ) 


PITCH  ANGLE 

S3-2  LAPE  HULA  1589  03/27/76  80  SECONDS  DATA 

11/17/78  PLOT  NO.  A  UT=19.12  LT=17.78  MT=22.30 

o  **.5  MEV  A  1.2  MEV  +  .2  MEV  X  .1  MEV  $  75.  MEV 

^  37.5  MEV  x 19  MEV  X  10.5  MEV  Y  6.8  MEV 


RTH 


PITCH  ANGLE 

S3-2  LAPE  GUAM  1592  03/27/76  120  SECONDS  DATA 

11/17/78  PLOT  NO.  H  UT=  .19  LT=19.W  MT=22.A1 

o  4.5  MEV  A  1,2  MEV  -f  .2  MEV  X  .1  MEV  $  75.  MEV 

A  37. 5  MEV  x  19  MEV  Z  10.5  MEV  Y  6.8  MEV 


■-i- > 
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5-2  LAPE  BOSS  1888  04/17/76  128  SECONDS  DATA 

11/16/78  PLOT  NO.  3  UT=  18.43  LT=  21.74  MT=22.86 

©  4.5  MEV  £  1.2  MEV  +  .2  MEV  X  .1  MEV  $  75.  MEV 

t  37.5  MEV  *19  MEV  *  10.5  MEV  Y  6.8  MEV 


NHR  T  m  ) 


PITCH  ANGLE 

13-2  LAPE  BOSS  1888  04/17/76  190  SECONDS  DATA 

11/16/78  PLOT  NO.  6  UT=18 . 51  LT=12.60  MT=6.72 

O  4.5  MEV  &  1.2  MEV  +  .2  MEV  X  .1  MEV  $  75.  MEV 

♦  37.5  MEV  *19  MEV  £  10.5  MEV  Y  6.8  MEV 


4 . 3  Evaluation  of  Generalized  Geometrical  Eactor 
by  Use  of  Monte  Carlo  Method 

In  the  course  of  the  analysis  of  the  p-a  and  LEPS  data 
a  reformulation  of  the  generalized  geometrical  Factor, 
(G(X,n))  for  these  instruments  was  performed  using  the 
Monte  Carlo  Technique. 

The  Generalized  Geometric  Factor  for  the  P-a  detector 
(in  the  coincidence  mode)  and  LEPS  was  found  to  be: 

(cf.  AFGL  Report  TR-78-0137) 


r 2"  j 

r»,  .r  ' 

G(n,A)  =  2r2  1 

i  1 

f  cost?  s i n‘?  j|  1-  (cos •  cost  -sir.',  sin‘3  sine) 

J 

■o 

ii 

o 

c 

'  'L 

‘  e-o 

(  1. 


1  ) 


where 


n  =  parametric  index  of  the  pitch-angle  distribution(sin  0  )°, 

X  =  angle  made  by  the  detector  telescope  axis  f  with  the 
geomagnetic  field  B,  or  spin-geomagnetic  angle, 

r  =  radius  of  the  detector  , 

0  =  polar  angle  made  by  the  particle  rath  vector  *p  and  X  , 
fjjj  -  Carl  -'j  (maximum  polar  angle  for  which  a  particular 
path  can  intersect  both  the  upper  and  lower  surfaces 
of  the  detector  telescope^, 

_  n/.imml.-il  .ingle  ol  the  enli-y  pom l  >'i  the  incident 
particle  in  the  detector  coordinate  system  as  shown 
in  Fig.  I  .  I, 

D  -  length  of  the  telescope  cylinder  . 

64 


_ 


_ _ - 


_ 


For  the  convenience  of  rofereice,  the  sketch  of  the  coordinate 
system  is  reproduced  in  Fig.  1.1. 

The  values  of  G(n,X)  as  a  fun  :tion  of  n  and  A.  have  been  evaluated 
bv  means  of  numerical  integration,  and  these  results  were  presented  in 
aforementioned  AFC. I  report  in  the  form  of  tables  and  graphs. 

. 

:  For  the  purpose  of  testing  the  results,  we  have 

re-evaluated  G ( n , A )  by  use  of  the  Monte  Carlo  Method1.  The  method  is 

1  2 

best  illustrated  by  applying  it  to  the  game  ler-ruiii  problem  . 

3 

It  has  been  known  that  a  number  of  approaches  have  been  devisee 

in  order  to  increase  the  accuracy  of  the  Monte-Carlo  evaluation  by 

(  reducing  the  variance,  or  the  standard  deviation,  of  the  calculations. 

f  1  x 

We  have  used  a  simple  example,  J  c'ds,  to  test  si:;  different  approaches. 

J0 

The  results  clearly  show  that  any  approach  can  be  chosen  so  lon.3  as  a 

I 

large  enough  sample  is  used.  Since  the  first  approach  ("mean  value" 
approach)  is  the  easiest  to  use,  we  chose  it  to  evaluate  C(n,A)  as  given 
by  (1.1).  The  re  ults  are  presented  in  this  report  and  compared  with 
the  results  of  AFG1,  Report  TR-78-0137. 
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4.3.1  Calculation  of  the  Generalized  Geometrical  Factor  G(n,A) 
by  Use  of  the  Monte  Carlo  Method  of  Integration 

In  order  to  apply  the  Monte  Carlo  Method  to  calculate  the  generalized 

geometric  factor  G(n,A)*,  we  first  need  to  rewrite  the  integral  formula, 

2 

Kq .  (1.1),  in  a  form  such  that  the  mean-value  approach  can  be  used  .  The 
reason  for  us  to  choose  this  approach  other  than  the  other  five  approaches 
to  the  Monte  Carlo  Method,  can  be  seen  in  AFGL-TR-78-0146 
which  shows  that  this  approach  is  the  easiest  to  use  and  its  accuracy 
can  be  brought  up  to  99%  when  several  thousand  random  numbers  are  used. 

The  suitable  form  deduced  from  Eq.  (1.1)  is  given  below: 

2  2t,6M  N 

C(n,A)  =  2r  — —  l  g  (W  (  2.1) 

i  =  1 


where 


S( 0 i » ^ i )  sinG^ 


_n/2 


1  -  (cosAcosG . -sinAsinG . sinA  . ) ' 

l  l  i 


cos 


tan0.]  tanO^  1  tanG.  J 


4' 


-J  > 


(  2.2) 


In  evaluating  G(n, A)  through  Eq.  (  2.1),  we  pick  two  random 

numbers,  say  p^  v^,  at  a  time.  Since  the  range  of  p^  and  v.  is  from 

0  to  1,  these  numbers  must  be  tran:  formed  to  0  and  A.  as  follows: 

1  1 


°i  =  Vi  (w 


here  ft  *  tan 
M 


1  2  r 


H 


(  2  .3) 


=  2nv^ 


The  computer  program  written  to  generate  G(n,X)  by  the  Monte  Carlo 
Method  follows.  The  results  are  given  in  Section  3. 
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4.3.2  Values  of  G(n,X)  as  Evaluated  by  Use  of  the  Monte  Carlo  Method 


For  a  comparison  of  the  values  of  G(n,X)  calculated  by  use  of  numerical 

integration  and  the  Monte  Carlo  Method,  we  list  the  results  for  n  =  5  for 

* 

various  X's  in  Table  3.1.  It  is  clear  that  even  with  3000  random  numbers 
(i.e.  N  =  3000),  both  results  agree  within  two  significant  figures. 

Excellent  agreement  is  reached  when  the  random  numbers  increase  to  100,000. 
See  Table  4.1  of  Report  TR- 78-0137  for  comparison. 

Here  we  indicate  two  aspects: 

(i)  For  a  simple  integration  in  one  or  two  dimensions,  the  results 
can  be  obtained  easily  by  either  method.  However,  for  a  complicated 
multi-dimensional  integral  for  which  the  numerical  integration  becomes 
increasingly  cumbersome,  the  Monte  Carlo  Method  might  turn  out  to  be 
the  only  workable  means  of  evaluation. 

(ii)  For  handling  data  of  a  statistical  nature,  the  use  of  the 
Monte  Carlo  Method  has  an  additional  advantage  because  it  also  yields 
the  standard  deviations. 


The  calculation  was  based  on  the  geometry  of  the  Narrow-angle  LEPS 
which  has  a  collimator  of  a  radius  of  0.347  cm  and  a  length  of 
3.175  cm.  Therefore  the  limits  of  0  are  from  0  to  tan  * (2r/D)  =  12.3°. 
The  limits  for  <J)  are  from  0  to  2n,  which,  are  based  on  the  assumption 
of  the  axisymmetry  of  the  particle  flux  with  respect  to  the  local 


geomagnetic  field. 
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Table  4.1  G(n,A)  Evaluated  from  Numerical  Integration  and  from 
Monte  Carlo  Method  for  n=5,  A  =  0°  to  90° 


C(5,A) 

(From  numerical 

A  integration  'j 

G(5,A) 

( From  Monte 
with 

0 

Carlo  Method 

N  =  3000  ) 

G(  5, A) 
(From  Monte 
with  N  *= 

0 

Carlo  Method 
100,000  J 

0 

. 4648E-6 

. 4815E-6 

.9212E-8 

. 4655E-6 

. 1580E-8 

10 

.  8903F.-5 

. 8931E-5 

. 2226E-6 

. 8906E-5 

. 3863E- 7 

20 

. 1014E-3 

. 9803E-4 

. 1831E-5 

.  1017E-3 

. 3238E-6 

30 

. 5218E-3 

. 5165E-3 

. 7616E-5 

. 5218E-3 

. 1315E-5 

40 

. 1655E-2 

. 1667E-2 

. 2043E-4 

.  1651E-2 

. 3535E-5 

50 

. 3789E-2 

. 3779E-2 

. 4153E-4 

. 3784E-2 

. 7282E-5 

60 

. 6818E-2 

. 6783E-2 

. 7010E-4 

. 6784E-2 

. 1219E-4 

70 

. 1011E-1 

. 10Q6E-1 

. 1Q13E-3 

. 1011E-1 

. 1728E-4 

80 

. 1269E-1 

.  1249E-1 

. 1221E-3 

.  1274E-1 

. 2108E-4 

90 

. 1367E-1 

. 1364E-1 

. 1292E-3 

. 1364E-1 

. 2265E-4 

Note  that: 


n  =  Parametric  index  in  the  theoretical  pitch-angle  distribution 
(sin  o'  ),n 

A  =  Angle  between  the  telescope  spin  axis  and  the  geomagnetic 
field  (Spin-geomagnetic  Angle), 

o  H  Standard  deviation  of  G(n,A), 
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4.4  Satellite  S72-1 


A  documentation  package  containing  a  tape  copy  of  our  S72-1 
data  base  in  the  form  of  a  latitude-longitude  sort  into  a  one  degree 
grid  preserving  the  pitch-angle  and  energy  spectrum  information  was 
prepared  for  distribution.  Specifically,  it  was  made  up  at  the 
request  of  Dr.  Vette  to  be  incorporated  into  his  next  version  of 
proton  model  AP4 ,  but  it  will  be  made  on  a  request  basis  to  qualified 
parties . 
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5.  e2-s2  high  energy  proton  spectrometer 

In  order  to  expand  our  capability  in  the  field  of  magnetospheric 
physics,  it  was  decided  to  study,  design  and  build  a  charged-particle 
(proton)  telescope.  A  literature  survey  was  begun  to  determine  the  latest 
techniques  used  in  the  construction  of  such  a  telescope.  We  also  initiated 
a  search  into  the  scientific  literature  to  determine  the  outstanding  problems 
concerning  the  trapped  proton  fluxes. 

It  was  found  that  a  proton  telescope  capable  of  measuring  small  proton 
fluxes  in  a  high  radiation  environment  of  high  energy  electrons,  in  the 
1-100  MeV  energy  range,  but  more  particularly  in  the  1-24  MeV  range,  should 
be  built.  Indeed,  the  completeness  and  reliability  of  the  existing  data 
covering  this  energy  range  for  protons  leave  much  to  be  desired.  This 
instrument  should  have  good  background  (high  energy  electrons,  bremsstrahlung) 
rejection  characteristics,  together  with  a  high  sensitivity. 

5 . 1  Background  Rejection 

The  background  problem  was  investigated  using  a  program  from  the 
Radiation  Shielding  Information  Center  at  Oak  Ridge  National  Laboratories. 

This  program,  called  "Tiger",  is  a  comprehensive  software  system  for 
solving  one-dimensional  multilayer  electron/photon  transport  problems. 

From  our  computer  simulations,  we  calculated  that  a  grated  shielding 
made  of  copper  and  lead,  would  be  very  adequate  for  the  measurements 
to  be  performed  in  the  heart  of  the  radiation  belts  where  substantial 
electron  fluxes  exist.  Indeed,  it  was  shown  that  when  an  electron 
spectrum,  as  shown  in  Fig.  5.1,  interacts  with  a  shield  of  adequate  thickness 
to  stop  100  MeV  Protons,  a  substantial  amount  of  bremsstrahlung  is 
generated  which  of  course  will  interact  with  the  Si (Li)  detectors  of  the 
telescope.  Such  simple  shielding  had  been  used  in  a  previous  proton 
telescope.  Figure  5.2  shows  the  normalized  transmitted  bremsstrahlung 


71 


spectrum  resulting  from  the  interaction  of  an  electron  input  spectrum 
as  shown  in  Fig.  5.1  (High  mode),  with  a  1.32  cm  copper  shielding 

2 

(Curve  1).  The  relative  intensity  is  normalized  to  1  electron/cm  -sec. 
5.3%  of  the  electrons  will  produce  bremsstrahlung  with  the  indicated 
photon  spectrum.  The  peak  of  the  distribution  occurs  at  ^300  keV. 

The  probability  of  interaction  of  these  photons  with  the  Si (Li) 
detector  is  rather  high  since  the  average  thickness  of  each  detector 
as  far  as  the  bremsstrahlung  is  concerned  is  several  times  greater 
than  the  actual  thickness  of  the  detector  because  of  the  average  angle 
of  incidence  of  the  photons.  The  second  curve  (2)  shows  the  transmitted 
bremsstrahlung  spectrum  through  the  same  1.32  cm  of  copper  followed  by 
1.20  cm  of  lead.  The  lead  is  added  to  absorb  preferably  the  photon 
peak  at  ^300  keV.  The  total  photon  spectrum  is  thus  reduced  by  a  factor 
of  ^32.  Furthermore  the  peak  energy  in  the  photon  spectrum  is  displaced 
to  ^900  keV  where  the  interaction  cross  section  with  silicon  is  reduced 
by  h  as  compared  to  the  300  keV  cross  section. 

At  low  altitudes  (few  thousand  of  km  or  lower),  where  the  electron 
flux  are  smaller,  but  high  energy  (E  >  100  keV)  proton  fluxes  are 
troublesome,  the  inclusion  of  an  active  shield,  together  with  the  passive 
shielding,  would  probably  be  advisable. 

Two  different  approaches,  as  far  as  the  detector  logic  configuration 
is  concerned,  are  being  considered.  The  first  one  makes  use  of  the  well- 
defined  energy-range  relationships  for  protons.  The  second  approach 
utilizes  the  relationship  between  incident  energy  and  energy  deposited 
in  a  given  thickness  of  silicon. 
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ENERGY  SPECTka 


ENERGY  OF  TRANSMITTED  PHOTONS  (MeV) 


S. 2  Range  -  Energy  Telescope 


5.2.1  Description  of  the  Instrument 

In  this  approach,  the  telescope  will  consist  of  n^  totally  depleted 
Si(Li)  detectors  of  various  thicknesses,  separated  from  one  another  by 
absorbers  of  appropriate  thicknesses.  Each  detector-absorber  pair 
determines  an  energy  channel  .  The  energy  of  an  incident  proton  will  be 
determined  by  the  last  detector  to  record  the  passage  of  that  particle. 

Figure  5.3  is  a  plot  of  the  range-energy  law  for  protons  in  aluminum. 

The  same  curves  can  be  used  for  silicon  since  the  ranges  of  protons  of 
equal  energies  in  the  two  materials,  different  by  little  more  than  1% 
between  6  and  100  MeV. 

As  one  example  of  this  type  of  telescope,  we  use  9  detector-absorber 
pairs,  thus  determining  9  energy  channels.  A  tenth  detector  provides 
the  anti-coincidence  to  eliminate  even  s  occurring  from  the  back  of  the 
detector,  as  well  as  an  additional  ene  gy  channel  100-300  MeV. 

Table  5.1  shows  the  different  energy  channel  characteristics  of 
the  telescope.  The  thickness  of  absorbers  following  the  detectors  is 
also  given.  Figure  5.4  shows  the  energy  deposited  vs  the  incoming  proton 
energy  relation. 

The  entrance  aperture  of  the  instrument  contains  a  sweeping  magnet 
which  will  deflect  electrons  with  E  —  1  MeV  away  from  the  silicon  detectors. 

The  telescope  will  utilize  totally  depleted,  transmission  type,  ring 
mounted,  silicon  surface-barrier  detectors. 

5.2.2  Descriptioi  of  the  Circuit  Logic 

Figure  5.5  is  a  block  diagram  of  the  circuitry  for  the  proposed  instrument. 
Each  detector  feeds  a  standard  modular  pre-amplifier,  amplifier,  and 
discriminator-shaper  in  cascade. 
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The  discriminators  in  channel  4  through  9  are  set  to  fire  at  300  keV, 
the  average  energy  deposited  by  a  100  MeV  proton  in  200  micron  depleted 
detectors;  the  thresholds  for  channels  1  and  2  are  set  at  60  keV,  and  that 
of  3  is  set  at  100  keV.  The  instrument  will  not,  therefore,  be  sensitive 
to  protons  above  300  MeV,  nor  to  electrons,  which  on  the  average  deposit 
about  14,  17,  35  and  70  keV  in  40,  50,  100  and  200^ m  of  silicon. 

The  shapers  in  channels  1  through  9  set  up  biases  within  the  diode 
steering  matrix  in  such  a  way  that  the  delayed  "write”  pulse  generated  from 
the  firing  of  the  channel  1  shaper  is  steered  only  to  the  one  output  wire 
corresponding  to  the  last  solid  state  detector  which  has  produced  a  pulse 
in  sequence.  Thus,  the  condition  for  a  signal  appearing  on  the  channel 
n  wire  is  an  n-fold  coincidence  of  pulses  in  all  detectors  1  through  n, 
and  a  (9-n)  fold  anti-coincidence  of  pulses  in  all  detectors  (n+1)  through  9. 
This  logic  eliminates  essentially  all  background  arising  from  particles 
(or  bremsstrahlung)  which  penetrate  the  shielding  on  the  sides  of  the 
detectors . 

Additional  AND  logic  is  provided  for  recording  protons  greater  than 
100  MeV  in  energy  (those  which  penetrate  all  detectors).  This  is  a  3-fold 
AND  gate,  requiring  signals  from  the  first  (1),  and  the  last  (9),  and  the 
central  (7)  detector  to  reduce  accidental  coincidences. 

The  outputs  of  the  steering  matrix  and  the  AND  gate  all  feed  counting 
logic.  In  channel  1,  the  counting  logic  is  a  standard  logarithmic  count- 
rate  meter  (LCRM)  and  a  x  16  Scaler  with  digital-to-analog  converter  (DAC) 
fed  by  a  single-sho*  circuit  (SS) 

The  LCRM  output  in  this  channel  will  allow  for  extremely  high  counting 
rates  which  might  otherwise  lead  to  ambiguous  interpretation  of  DAC  outputs. 
In  the  remaining  channels,  2  through  9,  each  counting  logic  consists  simply 
of  a  X32  scaler  with  DAC,  fed  by  a  SS.  The  scaler-DAC’s  are  mounted  in 
pairs  on  circuit  boards. 
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5.3  Energy  Deposited  -  Energy  Telescope 

5.3.1  Description  of  the  Instrument 

In  this  approach,  the  telescope  will  make  use  of  4  totally  depleted, 
transmission-type,  ring  mounted,  silicon  surface-barrier  detectors. 

Again,  a  fifth  detector  with  an  absorber  of  appropriate  thickness  in  front 
of  it  will  be  used  for  anti-coincidence  purposes  for  protons  with 
E  >  100  MeV,  and  events  occurring  from  the  back  of  the  telescope. 

Figure  5.6  shows  the  energy  -  incident  energy  in  such  a  telescope, 
consisting  of  4  detectors  of  thickness  30,  200,  1000  and  2000  pm. 
lable  5.2  indicates  the  energy  deposited  in  each  detector  for  given 
incident  energies.  The  energy  range  has  been  divided  into  10  energy 
channels 

5.3.2  Description  of  the  Circuit  Logic 

This  instrument  makes  use  of  similar  modules  described  in  the  previous 
section.  However,  now  the  output  voltages  from  each  detector  are  A/D 
converted.  These  voltages  are  then  compared  to  a  matrix  of  digitized 
voltage  interval  values  as  shown  in  Table  5.2.  When  following  an  event, 
the  proper  relationship  in  one,  two  or  three  detectors  is  shown  to  exist, 
the  proper  energy  channel  is  incremented. 
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Table 


Table  5.2 


Detector 


Energy  (MeV) 


Energy  Deposited 
(MeV) 


Energy  Deposited 
(MeV) 


Energy  Deposited 
(MeV) 


Energy  Deposited 
(MeV) 


6.  DETECTORS 


6.1  Detection  of  Protons  in  CR-39  Plastic  Track  Detector 

6.1.1  Abstract 

CR-39  plastic  is  used  as  a  detector  to  observe  monoenergetic 
protons.  Several  samples  of  CR-39  were  exposed  to  protons  of  energies 
1.5  MeV,  2.2  MeV,  3.2  MeV  and  4.3  MeV.  After  etching  measurements 
were  carried  out  on  track  diameters  produced  by  protons  in  all  samples, 
the  diameter  distributions  clearly  showed  an  excellent  energy  resolution 
of  protons  of  different  energies.  From  our  preliminary  analysis  of  the 
data,  it  appears  that  the  response  function  for  registering  protons  in 
CR-39  may  be  much  less  complicated  than  in  cellulose  nitrate. 

6.1.2  Introduction 

In  recent  years  Solid  State  Nuclear  Track  Detectors  have  found 
widespread  application  The  production  of  tracks  by  energetic  ions 

in  insulating  materials  is  a  widely-used  technique  for  detection  and 
identification  of  ions.  Their  utilization  showed  to  be  very  successful 
in  the  study  of  very  heavy  primary  cosmic  rays  and  the  recording  of 
fission  fragments.  Also,  there  have  been  some  investigations  exploring 
the  possibility  of  their  application  to  detect  protons.  Cellulose  Nitrate 
plastic  has  been  employed  as  a  detector  to  record  protons,  however,  CN 
suffers  from  being  inhomogeneous  and  anistropic  with  regard  to  its 
physical  characteristics.  These  defects  manifest  themselves  in  non- 
geometrical  track  profiles,  differences  in  sensitivity  in  a  given  sheet, 
and  varying  bulk  etch  rate.  The  use  of  plastic  sheets  cast  from  CR-39 
monomer  (allyl  diglycol  carbonate)  with  excellent  etching  properties, 
high  sensitivity,  and  high  uniformity  as  a  nuclear  track  was  reported 
recently  .  This  material  was  found  to  have  a  lower  detection  threshold 
(Z/0  =  9)  than  cellulose  nitrate  (Z/8  =  30)  and  a  smaller  variation  of 
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response  (M.%)  to  particles  of  a  given  ionization  rate  than  Lexan 
polycarbonate  (^3  to  8%).  Commercially  available  CR-39  is  capable  of 
recording  protons  of  1  MeV  and  below  as  well  as  6  MeV  alpha  particles. 

In  this  work  an  attempt  is  made  to  investigate  further  the  properties 
of  CR-39  as  a  nuclear  track  detector  with  particular  reference  to 
detection  of  protons. 

6.1.3  Experimental  Details 

Far  this  study  Polytech  CR-39  (1500  ym  thick)  was  used. 

Samples  of  size  (2. 5x2. 5  cm)  were  cut  out  and  exposed  to  a  beam  of 
protons  from  the  Van  de  Graaff  generator  at  Brookhaven  National 
Laboratory.  The  beam  was  tuned  to  four  different  energies,  viz., 

1.5  MeV,  2  2  MeV,  3.2  MeV  and  4.3  MeV.  The  corresponding  fluences 

5  2  6  2  6 

were:  2. 8x1.0  particles/cm  ,  0.44x10°  particles/cm  ,  0.67x10  particles/ 

2  6  9 

cm  and  0.91x10  particles/cm  .  All  irradiations  were  carried  out  in 
such  a  manner  that  particles  are  nearly  normally- incident  to  the 
surface  of  the  plastic  sheet.  After  irradiation  the  samples  were  etched 
in  a  solution  of  6.25N  Naoll  at  50°C  for  different  durations.  The  samples 
were  suspended  by  means  of  nichrome  wires  in  polyethylene  vessels 
containing  etchant.  Temperature  control  was  achieved  by  placing  the 
polyethylene  containers  into  a  regulated  water  bath.  Samples  were 
selected  from  all  four  irradiations  and  etched  in  four  different  batches 
for  7  hours,  17  hours,  30  hours,  and  48  hours.  In  order  to  distinguish 
the  background  against  tracks  of  protons,  always  a  virgin  sample  of 
CR-39  was  etched  with  etch  batch.  All  measurements  were  made  on  a 
Koristka  R4  microscope  using  a  80  Zeiss  objective  and  xlO  wide  field 
American  Optical  eyepiece.  A  total  of  1000  tracks  were  measured  to 
obtain  data  on  track  diameters. 
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6.1.4  Results 


In  Figures  1  and  2  are  shown  the  distributions  of  track  diameters 
for  protons.  The  sheets  were  etched  for  17  hours  (figure  1)  and  30  hours 
(figure  2).  In  the  case  of  17-hour  etch,  there  is  some  spread  in  track 
diameter  distribution  for  4.3  MeV  protons.  This  may  be  of  statistical 
nature  and/or  spread  in  the  beam  energy.  The  superior  properties  of 
CR-39  as  a  track  detector  are  clearly  demonstrated  in  figure  2.  Here  the 
energy  resolutions  for  protons  is  excellent.  Once  again  the  peak  in  track 
diameter  distribution  for  4.3  MeV  protons  is  very  broad.  Figure  3  shows 
etch  pit  diameter  as  a  function  of  proton  energy.  The  samples  etched  for 
48  hours  show  a  maximum  slope  thereby  indicating  a  better  resolution. 
Figure  4  shows  etch  pit  diameter  as  a  function  of  amount  of  bulk  material 
removed  from  one  surface.  From  our  preliminary  data,  it  appears  that  the 
dependence  of  etch  pit  diameter  on  the  amount  of  bulk  material  removed 
from  one  surface  seems  to  be  less  complicated  than  in  CN.  Particle 

identification  by  measurements  of  etch  pit  diameter  was  first  suggested 

31  41 

by  Somogyi  .  Recently  Somogyi  and  Szalay  discussed  the  kinetics 

of  track  diameter  growth  in  considerable  detail.  In  principle  the  method 

should  work  with  particles  incident  at  arbitrary  angles  on  a  solid 

surface,  but  in  practice  it  is  much  simpler  if  the  detector  can  be 

positioned  such  that  particles  are  nearly  normally  incident.  For  tracks 

with  large  cone  angles  such  as  those  of  protons,  the  diameter  is  a  more 

sensitive  function  of  ionization  rate  than  is  track  length.  Finally 

figures  5  and  6  show  tracks  of  1.5  MeV  and  2.2  MeV  protons  respectively. 

It  is  a  pleasure  to  thank  Drs.  Norman  Rohrig  and  Steve  Marino  of 
Brookhaven  National  Laboratory  for  help  with  proton  irradiations. 
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Figure  Captions 


The  diameter  distribution  of  the  etch  pits  of  protons 
(17  hour  etch) 

The  diameter  distribution  of  the  etch  pits  of  protons 
(30  hour  etch) 

Etch  pit  diameter  (D)  as  a  function  of  proton  energy  (E) 
The  parameter  on  each  curve  represents  etching  time. 

Etch  pit  diameter  (D)  as  a  function  of  amount  of  bulk 
material  removed  from  one  surface  (V^t) 

Tracks  of  1.5  MeV  protons  in  CR-39.  The  etching  time 
was  30  hours. 

Tracks  of  2.2  MeV  protons  in  CR-39.  The  etching  time 
was  30  hours. 
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7.  MAGAF  SYSTEM 

7.1  Systems  Generation 

Initial  preparations  have  been  made  to  do  another  systems  generation 
(sysgen)  on  the  Univac  minicomputer  which  services  the  AFGL  Magnetometer 
Network.  The  purpose  of  this  sysgen  is:  1)  to  allow  the  use  of  both  the 
teletype  and  the  Tektronix  CRT  at  the  same  time,  2)  to  incorporate  the 
DEC  printer  into  the  system,  3)  to  incorporate  the  new  drivers  into  the 
system,  i.e.,  the  magnetic  tape  and  DDU  drivers,  and  4)  to  incorporate 
the  existing  patches  contained  in  the  patch  file  into  the  system. 

To  accomplish  this  the  following  things  must  be  done:  1)  update 
the  current  systems  tape  to  reflect  the  new  drivers  and  to  include 
dummy  drivers  and  controller  tables  in  order  to  do  the  necessary  diag¬ 
nostic  testing,  2)  update  the  sysgen  directives,  and  3)  disassemble 
the  patch  file. 

The  two  new  drivers  have  been  modified  and  the  dummy  controller 
tables  and  drivers  have  been  made,  and  these  need  now  to  be  added  to  the 
system  tape.  The  controller  tables  for  the  DEC  printer  and  the  teletype 
have  also  been  written. 

7.2  Documentation  of  Minicomputer  Software 

Programs  written  for  use  on  the  Univac  minicomputer  which  services 
the  AFGL  magnetometer  network  are  in  the  process  of  documentation  and, 
in  some  cases,  overhauling.  Descriptions  and  use  of  those  documented 
to  date  are  given  in  this  section. 

For  all  minicomputer  programs  described,  operator  communication 
with  the  minicomputer,  is  through  the  Tektronix  CRT.  Similarly  all  out¬ 
put  from  the  minicomputer  is  received  on  the  CRT.  With  each  program 
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is  provided  the  commands  supplied  by  the  operator  and  the  output  from  the 
minicomputer,  both  of  which  appear  on  the  CRT  screen.  Hard  copy  of  the 
CRT  screen  Is  made  by  means  of  the  Tektronix  Hard  Copy  Unit,  which  is 
connected  to  the  CRT,  and  it  is  these  copies  which  are  used  to  illustrate 
the  programs. 

Some  programs  have  more  than  one  mode  of  operation.  Different 
options  can  be  selected  by  responding  at  the  CRT  to  prompts  made  by  the 
program.  Other  options  are  selected  by  means  of  the  sense  switches, 
referred  to  as  SSI,  SS2,  and  SS3,  respectively,  which  are  on  the  front 
panel  of  the  minicomputer.  Each  will  be  set  on,  off,  or  ignored  depending 
on  the  particular  program  and  the  option  desired.  At  other  times  the 
sign  switch  on  the  front  panel  can  be  used  to  signal  the  minicomputer. 
These  variations  are  explained  for  each  program  in  this  section. 

In  this  section  are  included  descriptions  of  the  following  programs 
and  procedures:  A)  TAPLOG,  B)  TAPDUP,  C)  Procedure  to  reconstitute  the 
system  from  tape  (DISKNU),  and  D)  NUFILE.  Examples  and  listings  are 
also  provided. 

7.2.1  Program  Descriptions  and  Procedures 

TAPLOC 

Program  TAPLOC.  reads  an  ARCIIIV  tape,  looks  at  the  time  of  the  first 
frame  of  each  record,  and  may  print  it  depending  on  criteria  explained 
below.  Some  diagnostics  arc  also  printed  in  the  case  of  defective 
records.  This  program  is  useful  for  finding  out  what  time  period  is  on 
a  particular  tape,  what  records  are  defective,  and  ,  in  conjunction  with 
other  programs,  for  finding  where  a  particular  time  period  on  a  tape 
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The  program  has  two  modes  of  operation,  the  difference  between 
them  being  the  amount  of  output  produced  by  the  minicomputer.  The 
modes  are  selected  by  means  of  sense  switches,  and  the  mode  can  be 
changed  while  the  program  is  running. 

With  SS3  on,  the  time  of  the  first  frame  of  every  record  is 
printed  whether  that  time  is  good  or  bad.  This  option  does,  however, 
produce  a  lot  of  output  since  there  are  usually  about  4000  records 
per  tape,  equivalent  to  4000  lines  of  output. 

With  SS2  on  and  SS3  off,  the  time;  of  the  first  frame  of  a  record 
is  printed  only  if  one  of  the  following  conditions  is  fulfilled: 

a)  The  time  of  the  first  record  is  always  printed. 

b)  The  time  is  printed  if  a  new  day  is  started,  i.e.,  the  time 
just  passed  midnight  on  that  record. 

c)  All  ends  of  file  (EOF's)  are  printed. 

d)  The  first  time  after  an  EOF  is  printed. 

e)  At  any  point  where  there  is  a  logical  gap  in  the  data,  the 
time  is  printed. 

f)  Had  times,  i.e.,  times  which  make  no  sense,  are  also  printed, 
e.g.,  the  day  is  greater  than  366,  the  time  contains  all  0's, 
etc. 

If,  while  running  with  only  SS2  on,  the  operator  desires  to  know 
what  the  frame  time  is  at  that  point  in  the  tape,  pressing  the  sign 
switch  on  the  front  panel  will  produce  the  current  time  being  processed. 
This  feature  is  handy  in  positioning  a  tape  at  a  particular  time  for 
another  program. 

Operator  commands  and  minicomputer  output  are  shown  in  Figure  7.1. 
The  operator  starts  the  program  by  typing  '/TAPL0G'.  The  program  asks 
which  magnetic  tape  drive  contains  the  input  tape,  "ENTER  INPUT  UNIT  : 
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1,2,3  FOR  Ml,  M2,  M3."  The  operator  responds  accordingly.  The  next 
line  reminds  the  operator  of  the  use  of  the  sense  and  sign  switches 
with  this  program. 

Following  this  the  times  are  printed  when  required.  The  format 
of  a  line  reading  across  is: 

a)  'SEQ'  means  sequence; 

b)  the  time  of  the  current  record, 

c)  'FOLLOWS', 

d)  the  time  of  the  previous  record, 

e)  the  number  of  the  previous  record  (or  the  number  of  the 
current  record  if  you  count  from  0). 

The  previous  record  is  printed  because  in  the  event  of  a  bad  time  one 

wants  to  know  where  on  the  tape  it  occurred.  The  format  of  the  times 

is  best  explained  by  example.  In  Figure  7.1,  the  time  of  the  first 

record  is  '8  94  15:25:40',  which  is  year  1978,  Julian  day  94,  15  hours, 

25  minutes,  and  40  seconds. 

Diagnostics  are  printed  in  the  case  of  defective  records.  'PAR' 
indicates  a  parity  error  on  that  record.  Note  that  this  could  be  due 
to  a  hardware  problem  as  well  as  a  bad  record. 

'BUM',  meaning  a  'bum'  record,  indicates  that  the  consistency 
check  made  on  the  record  failed,  i.e.,  the  program  could  make  no  sense 
of  the  record.  In  the  consistency  cheek  the  program  reads  through  the 
record  frame  by  frame,  reading  from  each  frame  the  length  of  that 
frame,  and  proceeding  to  the  next  frame.  If,  at  the  end  of  the  record 
the  count  does  not  make  sense,  the  record  is  declared  bad,  or  'BUM'. 

Note  that  TAPLOG  does  not  check  each  frame  time  within  a  record,  i.e., 
it  is  possible  to  have  bad  frame  times  within  a  record  which  TAPLOG 
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considers  as  'good'. 

Looking  at  Figure  7.1,  we  see  that  the  tape  started  at  the  time  just 
described  above.  Midnight  was  passed  after  record  #859  and  again  after 
record  #3203.  There  are  bad  times  on  records  #1824  and  #3162.  An  EOF 
follows  record  #2727.  The  rest  of  the  times  represent  small  gaps 
(logical  gaps)  in  the  data. 

A  listing  of  TAPLOG  is  shown  in  Figure  7.2. 

TAPDUP 

Program  TAPDUP  is  used  to  copy  selective  records  of  an  ARCHIV 
tape  on  to  another  tape.  It  is  usually  used  in  conjunction  with 
TAPLOG  since  it  distinguishes  records  on  the  input  tape  by  their  record 
number  on  that  tape  rather  than  by  their  frame  times.  However,  no 
TAPLOG  is  needed  if  one  wants,  for  example,  to  duplicate  the  whole  or 
a  part  of  a  tape.  Selective  portions  of  many  tapes  e.g.,  magnetic- 
storm  commencements,  can  thus  be  merged  on  to  one  tape. 

Operator  commands  and  minicomputer  output  for  TAPDUP  are  shown 
in  the  example  given  in  Figure  7.3.  The  operator  starts  the  program 
by  typing  '/TAPDUP'.  The  programs  asks  which  magnetic  tape  drive  con¬ 
tains  the  input  tape,  in  a  manner  like  TAPLOG.  When  the  operator  has 
responded,  the  program  asks  which  magnetic  tape  drive  contains  the 
output  tape.  After  the  operator  has  responded,  the  program  is  ready 
to  receive  commands,  whose  descriptions  and  formats  are  explained 
below: 

C  Copies  one  physical  record. 

S  Skips  one  physical  record. 

Cn  Copies  through  record  n. 
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C9999  Copies  to  the  end  of  tape. 

Sn  Skips  through  record  n. 

W  Writes  an  EOF  (unless  the  last  record  written  is  an  EOF, 

i.e.,  it  will  avoid  a  double  EOF). 

T  Terminates,  i.e.,  it  writes  a  double  EOF  and  exits. 

Q  Quits,  i.e.,  it  exits  letting  the  next  input  tape  be 

mounted.  This  is  used  when  parts  of  two  different  tapes 

are  to  be  merged  on  to  one  tape.  TAPDUP  must  then  be 

reinitiated  as  explained  above. 

Note  that  'n'  denotes  the  record  number  through  which  one  copies  or  skips 
(not  the  number  of  records) .  Note  that  the  record  number  corresponds  to 
the  record  numbers  provided  in  the  right-hand  column  of  TAPLOG,  where 
this  record  number  corresponds  to  the  frame  time  given  by  the  rightmost 
time,  i.e.,  the  time  of  the  previous  record.  For  example,  in  Figure  7.1 
record  number  859  starts  at  23:19:30. 

In  this  hypothetical  example  of  TAPDUP  shown  in  Figure  7.3  we  want  to 
merge  parts  of  two  different  tapes  to  make  one  tape.  Consider  Figure  1 

in  which  ARCHIV  tape  #8094  was  run  on  TAPLOG.  Suppose  we  want  to  com¬ 

bine  everything  after  the  first  EOF  (at  record  #2728)  with  the  first  50 
or  so  records  on  the  next  ARCHIV  tape,  #8096  whose  TAPLOG  is  not  shown. 

We  also  want  to  delete  the  bad  time  at  record  #3162.  Note  that  if,  in 
trying  to  copy  a  specified  period  of  time,  we  want,  to  be  more  precise 
as  to  what  frame  times  we  copy,  we  would  have  to  use  the  sign  switch  or 
SS3  to  find  out  exactly  at  what  records  the  beginning  and  end  times  are. 

An  explanation  of  the  procedure  for  this  example  is  outlined  below: 

S2728  Skip  through  the  EOF. 

C3161  Copy  up  to  the  bad  time. 

S  Skip  one  record,  i.e.,  the  bad  time. 
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C4007  Copy  up  to  the  EOF. 

Q  TAPDUP  exits;  no  EOF  is  written.  Control  returns  to 

JCP. 

/REW,M3  Rewind  unit  M3. 

Dismount  tape 

Mount  new  tape,  # 8096,  on  the  same  tape  drive. 

/TAPDUP  Re-enter  TAPDUP,  specifying  input  and  output  units. 

C50  Copy  the  first  50  records. 

T  TAPDUP  writes  a  double  EOF  and  exits. 

The  finished  tape  was  run  on  TAPLOG,  and  the  result  is  shown  in  Figure  7.4. 

A  listing  of  TAPDUP  is  shown  in  Figure  7.5. 

DISKNU 

DISKNU  is  a  restoration  of  the  system  from  tape. 

In  the  event  that  the  operating  system  has  to  be  put  back  on  the 
disc  from  magnetic  tape,  a  specified  procedure  must  be  followed.  This 
procedure  is  shown  step  by  step  in  Figure  7.6. 

Note  that  in  the  present  configuration  the  operating  system  resides 
on  the  bottom  platter  of  the  disc.  The  normal  boot  procedure  is,  with 
the  switch  above  the  front  panel  set  to  'DISC  BOOT',  to  press  successively 
'RESET',  'STEP/RUN',  and  'BOOT'  on  the  front  panel. 

NUFILE 

Program  NUFILE  ?s  a  utility  program  which  allows  the  user  to  go 
easily  from  magnetic  tape  to  disc  and  vice  versa.  It  dumps  out  or 
reads  in  partitions,  file  by  file,  thereby  allowing  selective  additions. 
Merging  and  garbage  collection  can  also  be  done. 

NUFILE  can  be  used  to  back  up  the  disc,  including  the  bottom 


platter  on  which  resides  the  system.  (See  the  system  recovery  procedure 
included  in  this  section.)  A  backup  copy  on  magnetic  tape  of  the  whole 
disc  is  now  routinely  kept. 

The  commands  which  NUFILE  provides  are  listed  below.  The  following 
nomenclature  is  used  to  specify  the  command  structure: 

a)  'unit'  or  'unit^'  is  always  a  disc-partition  logical-unit  name 
or  number.  (It  is  recommended  that  one  not  mix  logical-unit  names  and 
numbers . ) 

b)  'key'  is  the  protection  key,  if  any,  of  that  disc  partition. 

c)  'unit2'  actually  refers  to  the  tape  unit,  but  it  is  specified 

by  the  logical -uni*-  name  or  number  of  the  disc  partition  which  was  speci¬ 
fied  when  that  part  of  the  tape  was  made.  Note  that  to  go  from  tape  to 
disc  the  tape  must  have  been  previously  created  by  NUFILE. 

The  commands  are  as  follows: 

T  This  specifies  the  magnetic  tape  unit  for  dumping  or  partition 
loading.  The  default  is  unit  39,  DU. 

R  Rewind  the  tape  unit. 

D  D.unitj [,key[,unit2]]  dumps  unit^  to  unit2- 

L  L,unitj [ ,key [ ,unit2] ]  loads  unitj  from  unit2,  with  garbage  col¬ 
lection,  e.g.,  if  we  want  to  load  disc  partition  DO  from  that 
partition  named  D2  on  the  tape,  we  would  type  'L,D0,,D2'. 

G  This  does  garbage  collection  on  a  disc  partition;  the  magnetic 
tape  is  used  as  a  buffer.  'G.unit [ ,key] '  is  equivalent  to: 

D, unit [.key] 

R 

L.unit [ ,key] 

M  This  merges  files  on  unit2  on  to  unit},  if  they  are  not  already 
there.  The  format  is:  M.unit j [ .key [ ,unit2] ] . 

A  This  adds  files  from  unit2  to  unitj.  It  prompts  operation  for 
the  file  list.  Files  to  be  added  must  not  have  names  which 
correspond  to  files  already  on  unitj.  If  a  file  name  is  used 
which  conflicts,  the  program  will  print  'ALREADY  THERE'  and 
ignore  the  request.  The  format  is:  A, unit]  .key  ,unit2. 


Q  Quit,  i.e.,  exit  program  NUFILE.  Note  that  a  '/'  will  not 
bring  you  back  to  JCP;  NUFILE  considers  this  a  bad  request. 

NUFILE  can  be  used  to  retain  object  modules  across  systems  genera 
tions  (sysgens).  Note  that  in  the  sysgen  software  provided  by  Univac, 
each  module  must  be  run  through  LMGEN  every  time  a  new  sysgen  is  done. 
After  the  tape  unit  has  been  specified  etc.,  the  necessary  commands 
are : 

D ,  OM ,  D 
D,BL,E 
D,FL,F 

When  the  new  sysgen  is  done,  these  partitions  are  reloaded  through 
NUFILE  by  typing  'L,OM,D',  etc. 

The  program  listing  for  NUFILE  is  shown  in  Figure  7.7. 

7.3  Software  Written  for  the  CPC  6600 

Three  subroutines  relating  to  the  ARCHIV  tapes  have  been  written 
in  Fortran  for  use  on  the  CDC  6600.  They  were  written  to  be  used  with 
and  as  an  enhancement  for  the  existing  software.  The  listings  are  in¬ 
cluded  in  Figures  7.8,  7.9  and  7.10. 

The  first  routine  called  TIMECHK  will  check  the  time  and  station 
numbers  of  each  frame  and  will  print  appropriate  diagnostics  if  any 
time  is  incorrect  or  if  there  is  a  large  time  gap  between  frames. 

The  second  routine  called  TOTDAT  is  used  in  conjunction  with 
TIMECHK.  It  constructs  the  year,  Julian  day,  hour,  minute,  second, 
and  second  within  that  year  for  each  frame  from  the  array  IDATE  which 
is  provided  by  the  existing  software.  For  each  record  the  frame  time 
is  stored  in  a  series  of  arrays  called  YEAR,  JDAY,  etc.  by  an  index 
whose  value  corresponds  to  that  frame's  position  in  that  record.  For 


example,  for  frame  #4  in  a  given  record,  the  time  is  found  in  YEAR(4) , 
JDAY (4) ,  etc. 

The  third  routine  DUMP2T  takes  the  frame  data  which  is  in  a  245 
word  array  and  prints  it  in  a  format  which  is  the  same  as  that  produced 
by  program  DUMP2T  on  the  Univac  minicomputer.  The  purpose  of  writing 
this  routine  was  that  people  looking  at  the  data  would  not  be  required 
to  get  acquainted  with  another  format. 
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0001  DIMENSION  IA  ( 2561 , 2  ) 

0003  INTEGER  *BUFF,*UITCH 

0003  DIMENSION  ITIME ( S  ), JTIME ( 5  ) 

0004  IEOF-0 
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0037  IF<NL-NL1  >61,62,57 
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Fig.  7.2  (cont. 


/TAPDUP 

ENTER  INPUT  UNIT!  1,2,3  FOR 
“enter  OUTPUT  UNIT  1,2,3;  BUT  NOT  INPUT 


Fig. 


ENTER  INPUT  UNIT i  1,2,3  FOR  HI,  1*12,  R3 
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0035  INREC-0 

0036  DO  15  1-1,4 

0037  IQQ-ISUA8(IDIG(I  )  ) 
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TO  RESTORE  THE  BOTTOM  DISC  FROM  TAPE 


1.  Mount  the  magnetic  tape  labeled  "SAVED  SYSTEM  BOTTOM  DISC  5/4/79" 
on  logical  unit  DU. 

2.  Set  the  switch  above  the  CPU  to  "TTY  BOOT".  (It  is  normally  on  the 
"DISC  BOOT"  position.) 

3.  On  the  CPU  console,  press  RESET,  STEP/RUN,  and  BOOT.  The  system  will 
now  hang  waiting  for  the  teletype  to  come  ready. 

4.  Make  sure  that  the  teletype  is  connected  to  the  back  of  the  console, 
is  turned  on,  and  is  on  line.  Mount  the  paper  tape  labeled  "SYSTEM" 
RECOVERY  BOOT"  (about  18  inches  long)  so  that  the  first  frame  is 
exactly  at  the  read  point. 

5.  Turn  on  the  paper  tape  switch.  The  paper  tape  will  be  read  and  the 
Tektronix  will  beep.  A  small  part  of  the  magnetic  tape  will  be  read. 
The  CRT  display  will  show 

001032  (001000)  001000  000004  007020  000000 

DG**  DISKNU 

6.  On  the  Tektronix,  enter:  G1032.  RETURN.  The  program  will  halt  with 
the  program  counter  =  01513 

7.  Press  START  on  the  CPU.  About  1/2  to  1  inch  thickness  of  magnetic 
tape  will  then  be  read  in;  program  DISKNU  is  reading  in  from  tape. 

When  this  is  finished,  the  CPU  will  be  in  blinking  run. 

8.  Turn  the  switch  above  the  CPU  to  the  "DISK  BOOT”  position. 

9.  Boot  the  system  in  the  usual  way  by  pressing  RESET,  STEP/RUN,  and  BOOT. 
The  teletype  will  print  out"PATCHING  COMPLETED". 

10.  The  system  will  come  up  in  the  usual  way,  making  the  query: 

UNIT  1  OR  2?  (0  FOR  NO  ARCHIV) 

Resume  archiving  if  desired;  restoration  is  complete. 

11.  Turn  off  the  teletype  and  dismount  the  saved-system  tape. 


Fig.  7.6 
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0261  11  IF(I UHERE <2)~IDIR(4/ NT RY))12, 13/14 

0262  14  INEU-1 

0263  GO  TO  17 


0264  13  INEU-0 

0265  17  MH-NTRY 

0266  Pl-NTRY+l 

0267  16  IENTRY-IDIR(7,Mf1)  +  INEU 
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0337  3  CONTINUE 

0338  CALL  BORDIK IDIR,IDIR( 1,NTRY+1  ),NTRY  ) 

0339  CALL  0UTIT(IDIR,IDIR(1,NTRY+1),NFILES) 


0340  RETURN 

0341  END 
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8.  MICROPULSATION  ACTIVITY 


8.1  Micropulsation  Events 

The  reader  of  this  report  is  referred  to  AFGL-TR-78-0312  for  a  comprehensive 
treatment  of  Mathematical  Analysis  of  data  from  the  MAGAF  Network  including 
studies  of  Magnetic  Pulsation  Events  and  Micropulsation  Events. 

A  study  has  been  undertaken  to  determine  how  best  to  use  data  from 
the  MAGAF  Network  to  specify  and  predict  status  of  the  magnetosphere  by 
means  of  micropulsations. 

As  a  first  step  a  number  of  magnetically  quiet  days  on  which  all 
seven  stations  were  active  has  been  selected  for  study  to  establish  the 
quiet  time  behavior  of  micropulsations.  A  list  of  these  days  below, 
and  some  representative  magnetograms  for  28  May  78  (Figures  8. 1-8.8) 
are  given. 

9  April  78 

15  April  78 

16  April  78 

27  May  78 

28  May  78 

Data  have  been  taken  from  the  Varian  System  and  inputted  to  the 
CDC  6600  System  at  AFGL. 

Modifications  have  been  made  to  the  existing  program  for  filtering 
the  data;  this  filtering  is  now  in  progress. 

Maximum  entropy  program  will  be  used  to  obtain  the  frequency  spectrum 
of  the  filtered  data. 

After  obtaining  the  frequency  spectrum,  results  will  be  presented 
in  the  form  of  a  dynamic  spectrum. 

A  search  was  carried  out  of  the  network  data  on  the  5  quietest  days 
of  the  month  up  to  December  1977.  From  this  search  9  days  were  selected 
as  promising:  4  May,  27  Apr,  26  May,  27  July,  5  Sept,  9  Oct,  23  Nov, 
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24  Nov,  8  Dec  and  18  December.  However  these  proved  unsatisfactory  because 
of  gaps  in  the  data  or  noisy  data.  The  search  will  be  resumed  for  the 
period  when  all  seven  stations  are  ope  'ating. 

A  talk  was  given  at  the  spring  19  8  meeting  of  the  AGU.  The  abstract 
is  attached  below. 


LONGITUDINAL  EXTENT  OF  OCCURRENCE  OF  Pc4 

AND  Pc5  GEOMAGNETIC  PULSATIONS 

W.  Bellew  (Physics  Research  Division,  Immanuel 

1. 

2. 

023872  FOUGERE 

1978  Spring  Meeting 

College,  Boston  MA  02115) 

V.L.  Patel*  (Dept,  of  Physics,  Massachusetts 

3. 

Solar-Planetary  Relations 

Institute  of  Technology,  Cambridge  MA  02139) 

P.  Fougere  (Air  Force  Geophysics  Laboratory, 

4. 

Micropulsations 

Hanscom  AFB  MA  01731) 

5. 

No 

We  have  analyzed  magnetic  pulsation  events 
from  the  AFGL  magnetometer  network  located  at 

6. 

No 

^55°N  geomagnetic  latitude  and  covering  3  hours 
of  local  time  sector.  Pulsation  events  with 

7. 

None 

period  greater  than  80s  have  been  studied  on 

June  3,  1977.  The  results  indicate  that  the 
pulsations  of  Pc4  with  small  amplitudes  (few 
gammas)  are  localized  in  a  longitudinal 
sector  of  i  or  2  hours  extent.  However,  the 
pulsations  woth  large  amplitudes  (^10y)  and 

8. 

Bill  to: 

Emmanuel  College 

Physics  Research  Div. 

400  The  Fenway 

Boston  MA  02115 

longer  periods  are  observed  over  a  longitudinal 
extent  of  3  hours.  Detailed  results  of  wave 

9. 

1143 

characteristics  will  be  presented  and  will 
be  compared  with  theoretical  models  for  the 
pulsations . 


As  a  continuation  of  the  study  of  Micropulsation  Events,  the  program 
BINNER  was  written.  This  program  divides  the  frequency  range  over  which 
a  frequency  spectrum  has  been  taken  in  to  18  bins,  selects  the  5  largest 
amplitudes  and  orders  tnem  according  to  frequency.  If  one  of  these 
amplitudes  is  found  in  a  particular  bin,  its  order  is  printed  in  that 
bin.  For  instance  at  Washington,  (KST=1)  at  midnight  UT  (8148000000 
i.e.  day  148,  1978  at  0000  UT)  from  the  output  we  sec  that  the  largest 
amplitude  was  in  bin  5,  the  next  largest  in  bin  9  and  so  on.  If  none 


of  the  five  largest  ampliudes  was  found  in  a  particular  bin  then  a  zero 
is  entered.  After  the  location  of  the  amplitudes  in  their  bins  the 
power  in  the  spectrum  is  printed  out,  followed  by  the  magnitude  of 
the  top  5  amplitudes. 

8.2  High  Time  Resolution 

In  order  to  study  micropulsations  with  high  time  resolution  we  have 
undertaken  writing  a  plotting  routine  that  plots  the  search  coil  data. 
Using  such  plots  it  will  be  possible  to  pinpoint  the  time  at  which  a 
change  in  the  magnetic  field  occurs.  Furthermore,  since  the  search 
coils  measure  the  low  frequency  component  of  the  magnetic  field 
is  effectively  filtered  out.  This  feature  is  illustrated  by  a  comparison 
of  the  following  figures  (Figures  8.9-8.14)  which  show  24  hours  of  search 
coil  and  fluxgate  data  for  the  same  periods. 

This  program  is  flexible  enough  so  that  it  will  be  able  to  plot 
for  a  period  as  small  as  1  minute. 
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8.3  Sudden  Storm  Commencement  (SSC) 

A  talk  was  given  at  the  AFGL  Magnetometer  Network  Workshop  on  the 
AFGL  network  data.  Plots  of  two  SSC  events  were  presented  to  give  an 
idea  of  what  the  data  are  like  as  recorded  by  the  fluxgate  and  search 
coil  magnetometers.  It  was  pointed  out  that  the  SSC  stands  out  more 
clearly  in  the  coil  data  than  in  the  fluxgate  data.  Comparisons  were 
made  in  all  three  components.  The  data  showed  on  of  the  SSC's, 

10  April  1978,  propagating  from  east  to  west,  away  from  the  subsolar 
point. 

Eighteen  SSC's  in  1978  were  selected  from  the  network  data  for 
study.  Approximately  one  hour  of  data,  starting  about  15  minutes 
before  the  SSC,  were  copied  onto  a  scratch  tape  from  the  original 
archive  tapes,  and  from  there  copied  onto  a  disk  pack  for  easier 
accessibility.  A  list  of  the  starting  times  of  the  SSC's,  given 
in  the  format,  YDAYHRMISC,  is 

8025152020 

8056191630 

8067142530 

8100125020 

8107233000 

8120094300 

8130195520 

8141023010 

8149182030 

8152213020 

8153090010 

8180093950 

8180204930 

8192103940 

8193000120 

8194000040 

8194231510 

8352002120 


Plots  of  an  SSC’s  trace  are  shown  for  8  March  73  at  five  stations 
(Figures  8.15-8.19).  The  data  plotted  here  are  from  the  search  coils. 
Since  the  search  coils  measure  changes  in  the  field  every  fifth  of  a 
second,  it  will  be  possible  to  study  these  events  with  high  time 
resolution. 
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